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SUMMARY: Pressure measurements were made at one section of an aerofoil 
of constant chord spanning a wind tunnel. A flat plate was mounted in the 
tunnel upstream of the aerofoil, so that the wake from the plate produced a 
spanwise velocity gradient at the aerofoil. By traversing the plate across the 
tunnel the position of the wake could be altered relative to the pressure holes 
in the aerofoil, so that the complete pressure distribution over the aerofoil 
in the neighbourhood of the wake could be found. The results showed that, 
because of the secondary flow produced by the interaction between the aero- 
foil and the wake, the variation of pressure across the wake for a given 
chordwise position was very small, except close to the leading edge. Thus 
the lift coefficient, based on the velocity outside the wake, was nearly constant 
across the wake. The form drag coefficient, based on the same velocity, was 
considerably reduced in the wake. Even when the form drag coefficient was 
based on the local velocity of the stream, a considerable reduction was 
found in the wake. 


1. Introduction 


The experiments described in this paper were made to investigate the flow past 
an aerofoil, in a stream whose velocity varies in the spanwise direction. Such a 
flow occurs, for example, at the wing root of an aircraft, where a part of the wing is 
within the boundary layer on the fuselage, and in experiments on aerofoils spanning 
wind tunnels, where the flow at the ends of the wing is affected by the boundary 
layers on the tunnel walls. 


There have been several attempts to investigate this problem theoretically, 
neglecting the effects of viscosity. In the earlier work, including that of Karman and 
Tsien", the assumption was made that the velocity gradient in the spanwise direction 
was small compared with the chordwise velocity gradient over the aerofoil, so that 
the flow at each section of the aerofoil could be considered as two-dimensional. In 
the experiments described here the spanwise velocity gradient is comparatively large, 
so that in any theoretical analysis of the problem the foregoing simplifying assump- 
tion cannot be made. 


In a theory developed by Hawthorne’, a small perturbation velocity is super- 
posed on a basic quasi two-dimensional flow to form the complete three-dimensional 
flow. The assumption of a small perturbation velocity means that the theory cannot 
be applied to an aerofoil with a rounded leading edge, having a stagnation point. 
(Hawthorne has used the theory to calculate the flow past an aerofoil with a cusped 
leading edge, at zero incidence.) 
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Mendelsohn and Polhamus®) have measured pressure distributions on a wing 
protruding from the wall of a wind tunnel, at various distances from the wall and at 
two angles of incidence. At the higher incidence the results were complicated by 
non-uniform stalling, but at the lower incidence a comparatively small change of 
chordwise pressure distribution was found as the wall was approached. The lift 
coefficient, based on the velocity outside the boundary layer, was only reduced by 
about 10 per cent. at a station very close to the wall. The drag of the aerofoil was 
not investigated. 


SIDE OF TUNNEL 
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SIDE OF TUNNEL 
Fig. 1. 
Plan of aerofoil and plate in tunnel. 


The present experiments were made in order to extend the work of Mendelsohn 
and Polhamus. To simplify the experiments a wake, produced by a flat plate 
upstream of the aerofoil, was used instead of the tunnel-wall boundary layer. This 
enabled the aerofoil to be kept fixed in the tunnel, the wake being traversed along 
the span of the aerofoil by moving the upstream plate. The aerofoil was made with 
sufficient pressure holes at one chordwise station to enable the drag due to normal 
pressures to be determined, in addition to the lift. 


Some measurements on an aerofoil in a wake have also been made by Ling 
and the results have been reproduced by Hawthorne’. Ling’s experiments were 
confined to measurements of drag on a balance, however, and no measurements of 
pressure distribution were made. The differences of drag observed by Ling were 
very small, so that any small errors in the measurements would have rather large 
effects on the conclusions. 


NOTATION 


chord of aerofoil 


Cy form drag coefficient, based on velocity outside wake 

Cy, lift coefficient, based on velocity outside wake 

C, pressure coefficient, based on velocity outside wake 

Cy’ form drag coefficient, based on local velocity at half-chord position 
C,’ _ lift coefficient, based on local velocity at half-chord position 


AC,’ _ difference between values of Cy’ outside wake and at specified position 
within wake 
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difference between values of C, on lower and upper surfaces of aero- 
foil, at specified position on chord 


local velocity in wake in absence of aerofoil 

velocity outside wake in absence of aerofoil 

distance along chord from leading edge of aerofoil 

distance from chord line of aerofoil, measured perpendicular to chord 
value of (y/c) at point (near leading edge) where C,=0 

angle of incidence of aerofoil 


change of (y/c), associated with increase of incidence required to 
produce (in a uniform stream) observed increase of Cy,’ 


difference between values of (y/c), outside wake and at specified 
position within wake 


value of at 


2. Experimental Details 


The experiments were made in an open-return wind tunnel, with a working 
section 74 in. wide by 20 in. high, at a wind speed of about 120 ft./sec. The aerofoil 
chord was 5°5 in., giving a Reynolds number of about 3:5 x 10°. Fig. 1 shows a 
plan view of the aerofoil in the tunnel. Following a suggestion of Mr. W. D. 
Armstrong, a perforated flat steel plate was used to produce the wake; this provided 
a simple method of making a plate with a rough surface that would produce a 
relatively thick wake. The plate was 0-07 in. thick, 5 in. long in the direction of 
the stream, and 13-8 in. high. The perforations were of 0-097 in. diameter, spaced 
about 0-14 in. apart. 


Fig. 2. 


Section of N.A.C.A. 0015-64 aerofoil, showing positions of pressure holes. 


TABLE I 
ORDINATES OF N.A.C.A. 0015-64 AEROFOIL SECTION 
x/c yle x/c yle x/c yle 
0-0125 0:0227 0-15 0:0607 0:60 0:0664 
0-025 0:0307 0:20 0:0660 0-70 0:0560 
0-05 0:0407 0°30 0:0727 0-80 0°0415 
0-075 0:0476 0:40 0:0750 0:90 0-0233 
0:10 0:0530 0°50 0:0727 0-95 0-0127 


Leading edge radius=0-0247c 
Trailing edge thickness=0-0030 c. 
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Fig. 3. 
Velocity profiles in wake at position of aerofoil. 


The aerofoil was of N.A.C.A. 0015-64 section (symmetrical) and was fitted with 
30 pressure holes. The aerofoil section and the positions of the pressure holes are 
shown in Fig. 2 and the ordinates are given in Table I. The pressure holes were 
closely spaced near the leading edge and one was also fitted at the trailing edge, so 
that the lift and form drag coefficients could be determined accurately from the 
measured pressure distributions. 


The wake produced by the plate was traversed with a pitot tube in the absence 
of the aerofoil at positions corresponding to the leading and trailing edges and at 
the half-chord position. The resulting velocity profiles are shown in Fig. 3. In all 
calculations involving the velocity in the wake at the position of the aerofoil, the 
velocity at the half-chord position has been used, except in considering the pressure 
at the front stagnation point, where the velocity at the position of the leading edge 
has been used. 


3. Pressure Distributions 

Chordwise pressure distributions were measured at incidences of 0°, 2°, 5° 
and 8°, on the wake centre line and at distances of 4, 4, 1 and 14 in. on either side of 
the centre line. For each of these distances from the wake centre line, mean values 
for the two sides of the wake were used in plotting the pressure distributions. Some 
further measurements at an incidence of 10:5° were rejected, because the pressure 
distributions showed that outside the wake there was a laminar separation near the 
leading edge on the upper surface, whereas in the more turbulent flow inside 
the wake this separation was suppressed. No separation of this kind near the leading 
edge was found at any of the lower incidences, but in some conditions there was 
evidence of a laminar separation on the lower surface near the trailing edge. This 
had only a small effect on the pressure distribution and is discussed later in 
considering the lift. 
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Fig. 4(5). 
Pressure distributions at ,=8° plotted against y/c. 
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Fig. 5. 
Maximum pressure coefficient on aerofoil. 


For a given chordwise position on the aerofoil, it was found at all incidences 
that the pressure coefficient C,, based on the velocity outside the wake, showed little 
variation across the wake, except at positions close to the leading edge of the 
aerofoil. This is shown, for example, in Figs. 4(a) and (b), where pressure distribu- 
tions are plotted for <=8° at the centre of the wake and outside the wake. Fig. 4(a) 
shows that the small variation of pressure distribution near the leading edge can 
have only a very small effect on the lift per unit span. 


Figure 4(b), where C, is plotted against y/c, shows more clearly the change of 
pressure distribution near the leading edge due to the presence of the wake. It can 
be seen that there are two main effects; in the wake the pressure at the stagnation 
point is reduced, as would be expected, and also the stagnation point is displaced 
along the lower surface away from the leading edge. Similar results were obtained 
at lower incidences, but of course at zero incidence there could be no displacement 
of the stagnation point. 


For a given position in the wake the maximum value of C, on the aerofoil was 
found to be independent of incidence between 0° and 8°. In Fig. 5 this maximum 
value (C, max) is plotted against distance from the centre of the wake. The dotted 
curve shows (u/U)? in the absence of the aerofoil at the position of the leading edge. 
These curves show that C, max is slightly less than (u/U)? over the whole of the wake, 
suggesting that the secondary flow in the spanwise direction prevents the formation 
of a true stagnation point. The difference between the two curves at the centre of 
the wake is surprising, however, because by symmetry there cannot be any span- 
wise flow at the centre. 


The pressure distribution curves (as in Figs. 4(a) and (b)) were integrated by 
means of a planimeter to derive coefficients of lift and form drag. The maximum 
possible errors in C; and Cy were estimated to be about +0-005 and +0-001 
respectively. 
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Fig. 6. 
Variation of lift coefficient across wake. 
4. Lift 


Figure 6 shows the lift coefficient, based on velocity outside the wake, plotted 
against distance from the centre of the wake. At the centre of the wake the lift 
coefficient is only reduced by about 2 per cent., although the local velocity (at the 
half-chord position) is about 15 per cent. less than the velocity outside the wake. 
The change of lift coefficient is as small as this because the variation of pressure 
across the wake is confined to a small region near the leading edge, and this has only 
a very small effect on the lift coefficient (Fig. 4(a)). 


Figure 7 shows the difference of pressure coefficient between the lower and 
upper surfaces, for an incidence of 8°. This provides a further illustration of the 
small spanwise variation of lift across the wake, and shows that the small variation 
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Fig. 7. Fig. 8. 
Chordwise load distribution at ¢=8°. Chordwise load distribution at ¢=5°. 


of pressure across the wake cannot be considered to be caused by a change of 
effective wind speed. 


Figure 8 shows the chordwise load distribution at an incidence of 5°. A 
noticeable feature of this diagram is that outside the wake AC, becomes negative 
near the trailing edge, whereas this does not occur at the centre of the wake. It is 
believed that the negative values of AC, are associated with laminar separation on 
the lower surface at the rear part of the aerofoil, and that this separation is 
suppressed in the more turbulent flow in the wake. Similar effects have been 
observed previously on aerofoils with rather large trailing edge angles; for example, 
Preston, Sweeting and Cox'*) found negative values of AC, near the trailing edge 
of a Piercy 1240 aerofoil, which were eliminated when transition wires were fitted 
at O-lc. 


Figure 8 shows that, if the laminar separation just discussed had not been 
present, the variation of C, across the wake would have been slightly greater than 
that shown in Fig. 6. 


5. Drag 


Figure 9 shows the form drag coefficient, plotted against distance across the 
wake; the full lines refer to the coefficient Cp (based on velocity outside the wake), 
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Fig. 9. 
Variation of form drag coefficient across wake. 


while the dotted lines refer to Cy’ (based on local velocity in the wake). The values 
of Cp and C,’ are fairly high, because of the low Reynolds number. It is clear that 
both Cp and Cy)’ are considerably reduced in the wake, especially at the higher 
incidences. (If there were no secondary flow, and no induced upwash or downwash 
due to spanwise variation of circulation, Cp’ would be constant across the wake.) 


6. Discussion 


The flow considered in these experiments is a complicated one, but some of 
the main features can be explained qualitatively, as follows. 

Near the leading edge of the aerofoil the pressure is reduced inside the wake, 
because the total pressure of the stream is reduced, so that there is a spanwise 
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TABLE II 


DISPLACEMENT OF POSITION ON AEROFOIL AT WHICH C,=0 y 


Distance from 


centre of wake 0 0:25 0°5 
(in.) 
A 0:0024 0-0013 0-0003 
am2° 8 0:0030 0-0021 0-0008 
0-0017 0-0013 0-0005 
A 0-0041 0-0026 0:0008 
ams? 8 0-0035 0:0017 0-0003 
0:0022 0-0009 0-0000 
A 0:0053 0-0036 0-0012 
a=8° 4 8 0-0067 0-0043 0-0013 
0-0054 0-0035 0-0010 


pressure gradient causing a flow towards the centre of the wake. Slightly farther 
back on the aerofoil there are pressure gradients tending to cause a spanwise flow 
away from the centre of the wake, and in the resulting flow there is a component of 
vorticity in the direction of the undisturbed stream. The effect of the spanwise flow 
on the pressure distribution is to equalise the pressures along the span of the 
aerofoil, for a given chordwise position. Fig. 4(a) shows that this equalising process 
is very effective except at points close to the leading edge. 


Since the pressure distribution is almost unaffected by the wake except near the 
leading edge, there cannot be much variation of lift coefficient across the wake. 
There is, however, a considerable reduction of form drag, because the reduction of 
pressure near the leading edge has a significant effect on the chordwise force 
component (Fig. 4(5) ). 


It has already been mentioned that in the wake the stagnation point is 
displaced along the lower surface away from the leading edge (Fig. 4(b)). Ina 
uniform stream such a displacement of the stagnation point would be associated 
with an increase of circulation, and hence of lift coefficient, and it is of interest to 
compare the increase of lift coefficient, calculated from the displacement of the 
stagnation point, with the observed increases of the coefficient C,’, based on local 
velocity in the wake. 


From diagrams of the kind shown in Fig. 4(b) the exact position of the stagna- 
tion point cannot be found, but it is possible to determine fairly accurately the 
value of y/c at the point on the aerofoil (near the leading edge) at which C,=0. 
Denoting this value by (y/c),, let 5 be the difference between the values of (y/c), 
within the wake and outside it. Within the wake the lift coefficient C,’, based on 
local velocity, is increased. From the measurements made outside the wake it is 
possible to find the change of (y/c), (say A) associated with the increase of incidence 
required to produce this increase of local lift coefficient. The quantities 5 and A 
may then be compared. 
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TABLE IiIl 


NEGATIVE INDUCED DRAG 


Distance from 


centre of wake 0 0:25 0:5 
(in.) 
Cr’e 0 0 0 
0-0044 00021 0:0008 
00026 00013 0:0003 
ACp’ 00027 00020 0-0012 
0:0156 0:0091 0:0025 
ACp’ 0:0057 00041 00010 
0:0375 0-0237 0:0076 
ACp’ 0:0158 0:0113 0:0018 


The analysis is complicated, however, by the fact that 5 is not zero even at 
zero incidence, because of the effect of the wake on the pressure distribution near 
the leading edge. If 5, is the value of 5 at zero incidence, for a given position in 
the wake, it may perhaps be more appropriate to compare A with (8—6,). 


Table II shows values of A, 5 and (8—8,), calculated as just explained. The 
values given in the table show that A agrees in order of magnitude with either 5 or 
(8—5,), suggesting that there is, in fact, an increase of circulation round the aerofoil 
in the wake, corresponding roughly with the observed increase of C,’. This con- 
clusion is consistent with observations of wing stalling made in other wind tunnels, 
where it has often been found that the part of the wing inside the tunnel boundary 
layer stalls at a lower incidence than the remainder of the wing. 


The foregoing discussion of the variation of C,’ across the wake suggests that 
the observed reduction of form drag coefficient in the wake may be partly caused 
by an induced upwash due to the secondary flow, giving a “ negative induced drag.” 
Basing the coefficients of both lift and form drag on local velocity, the reduction of 
form drag coefficient due to this effect should be C,’¢, where ¢ is the angle of induced 
upwash in radians. Values of C,’< have been calculated from the measured lift 
distributions, assuming that the observed increase of C,’ in the wake is due entirely 
to the induced upwash «, and that there is no variation of lift curve slope (based on 
local velocity) across the wake. In Table III these values are compared with the 
observed reduction of form drag coefficient AC)’. At high incidences the observed 
reduction of form drag coefficient is considerably less than the calculated “ negative 
induced drag,” i.e. there must be some other effect giving a substantial increase of 
drag. At zero incidence the negative induced drag is necessarily zero, although a 
reduction of C,’ is still observed. Thus in this case there must be some other effect 
causing a reduction of drag, independent of the negative induced drag. 


Since form drag is entirely due to thickening or separation of the boundary 
layer at the rear of the aerofoil, it might be expected that a reduction of form drag 
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would result from the better mixing in the boundary layer, caused by the secondary 
streamwise vorticity, the mechanism being similar to that of vortex generators. 


The suggestions just made, to account for the observed reduction of Cy’, can 
only be regarded as tentative ones. It is clear that a complete explanation of the 
observed drag changes would be very complicated. Nevertheless, the fact that 
the observed reduction of Cy’ in the wake is greater at the higher incidences suggests 
that negative induced drag may play some part. 


7. Conclusions 


It has been shown that, when an aerofoil of constant chord is placed in a stream 
with a spanwise velocity gradient, there is only a small spanwise variation of pressure 
coefficient for a given chordwise position. Basing coefficients of lift and drag on a 
constant velocity, it is found that there is little spanwise variation of lift coefficient, 
but there is a considerable reduction of form drag coefficient in a region of reduced 
stream velocity. The observed spanwise variation in the position of the stagnation 
point is roughly consistent with the variation of lift coefficient based on local 
velocity. 


The foregoing results have some interesting practical applications. For 
example, when measurements are made on an aerofoil spanning the working section 
of a wind tunnel, the lift per unit span decreases only a little where the aerofoil 
enters the boundary layers on the walls of the tunnel. (This has already been shown 
by the work of Mendelsohn and Polhamus™’). Again, at the junction between the 
wing and fuselage of an aeroplane, there should be a reduction of local drag 
coefficient of the wing due to the boundary layer on the fuselage. In practice a net 
reduction of drag has never been achieved, probably because there is a counter- 
acting increase of drag due to the effect of the wing and its boundary layer on the 
fuselage boundary layer, but the results given here do at least suggest that the ideal 
to be aimed at in such cases is a reduction of drag. 
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Experiments with an Electrical Analogue for 
the Extension and Flexure of Flat Plates 


P. J. PALMER, M.Sc., Ph.D., D.I.C., A.R.C.S., A.F.R.Ae.S., A.Inst.P. 
and 


S. C. REDSHAW, D.Sc., Ph.D., M.LC.E., F.R.Ae.S. 


(Department of Civil Engineering, University of Birmingham) 


Summary: Details of the design, construction and use of an electrical 
resistance network for solving certain of the problems encountered in the 
flexure and extension of flat plates are given. Some typical experiments have 
been carried out with the aid of this network and the method of operation of 
the analogy is described. The results of these experiments are presented and 
they are compared with the corresponding theoretical values. The comparison 
is as good as could be expected in view of the limitations imposed by the very 
coarse mesh which was used; essentially the same results would have been 
obtained by the exercise of a relaxation computation, using the same size of 
net. No practical difficulties would be encountered in building a large 
network having a very fine mesh. 


1. Introduction 


An earlier report") considered the feasibility of constructing various electrical 
networks for the solution of the partial differential equations concerning the 
extension and flexure of flat plates. The most promising of these networks was 
considered to be a pure resistance network, and it was recommended in the report 
that such an analogue should be constructed and tested. 


The purpose of the present report is to give details of the construction of a small 
pure resistance network, to show how it can be used to solve the relevant differential 
equations and to give the results obtained with it. 


It must be emphasised that the network constructed and reported here was only 
a small one, so that all the results are necessarily those obtained by a solution 
appropriate to a coarse grid. If the results obtained with this small network are 
regarded as indicating the analogy to be satisfactory, then a larger network can be 
built without any difficulty. 


Received September 1954. 
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distance between mesh points 
flexural rigidity of plate 
electrical potential in high potential grid 
a continuous function 
neighbouring ordinates in the finite difference representation of f 
resistances in the gradient measurement circuit (Fig. 7) 
(M.+M,)/{D(1+y)} = —V?w 
bending moments in plate in x and y directions 
direction normal to plate boundary; integer 
intensity of distributed load normal to plate 
value of resistances connecting low potential grid to medium potential 
grid 
value of resistances connecting medium potential grid to high potential 
grid 
value of resistances in low potential grid 
electrical potential in low potential grid 
electrical potential in medium potential grid 
normal deflection of plate 
axes in plate 
stress function 
Poisson’s ratio 
value of resistances in medium potential grid 
potentials in the gradient measurement circuit (Fig. 7) 
stress components in plate 
3? 0? 
ax * 


ax * axtay 


position in grids corresponding to a typical point in the plate 


neighbouring points to that defined by suffix 0 (Fig. 1) 
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ELECTRICAL ANALOGUE 


2. The Mathematical Background 
2.1. VARIOUS CASES 


Extension and flexure problems of flat elastic plates are of very great 
importance in applied elasticity and, as shown by Southwell’, there is a definite 
interconnection between the two classes of problems. Details of the mathematical 
analogies which exist between the flexural and extensional problems have been 
admirably presented by Fox and Southwell’, who considered the problem of a flat 
elastic plate which is either bent or stretched. In all, four cases were presented 
since the edge conditions may specify either traction or displacements, both in the 
flexural and in the extensional problems. : 


The four cases may be defined as follows: — 


(A) The flexural problem in which displacements are specified at the 
boundary. 

(B) The flexural problem in which forces are specified at the boundary. 

(C) The extensional problem in which displacements are specified at the 
boundary and the only body forces are gravitational. 


(D) The extensional problem in which forces are specified at the boundary 
and the only body forces are gravitational. 


Normally Cases A and D are solved by one mathematical method and Cases B 
and C by another, although Southwell has shown that all four cases are, in fact, 
mathematically identical. To attempt to solve all four cases by means of one 
solution, either mathematically or electrically, leads in some cases to boundary 
conditions which are intractable. 


It is for this reason that the pure resistance network was used to solve problems 
of Cases A and D alone. Case A entails the solution of the equation 


V4w — q 
w= . . . . (1) 


such that usually w and its first normal derivative take specified values at the 
boundary; w is the normal deflection of the plate, g the intensity of the distributed 
load, and D the flexural rigidity of the plate. Case D entails the solution of the 
equation 


such that usually ¢ and its first normal derivative take specified values at the 
boundary, ¢ being a stress function. 


If derivatives higher than the first are specified at the boundary then it may still 
be possible to obtain a solution by means of the network, and further comments 


will be made about this later. 
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2.2. PARTICULAR FLEXURAL CASE 


Equation (1), which was first derived and considered by Lagrange, may be 
written as 


0? 0? \(0?w q 


Now the bending moments, in a flat plate, when expressed in terms of the 
normal deflection have the form 


0?w (4) 
aw —(M.+M,) _ 
Hence ax? +i =—M ‘ ‘ (5) 


the last equality being by definition. 


Thus the solution of equation (3) reduces to the solution of the two interrelated 
second-order partial differential equations given by 


0x? dy? ~D 
2 2 

= 0?’w _M 


2 = 


ox 


(6) 


and 


Any solution using an electrical resistance network is really analogous to the 
corresponding finite difference solution and not to the differential solution, although 
the difference between the two theoretical solutions can be as small as desired 
provided that the chosen ordinates are sufficiently close together. However, from 
the point of view of comparing the electrical network with theory, the finite difference 
forms of equation (6) are required. If the suffix 0 represents a typical ordinate 
defined by some value of x and y, and if the suffixes 1, 2, 3 and 4 represent neigh- 
bouring points on a square grid of side a, then the finite difference form of 
equations (6) becomes 


=0 


D 
and 


M,+M,+M,+M,-4M,+ 


the expression g, representing the intensity of loading at the chosen position. 
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2.3. PARTICULAR EXTENSIONAL CASE 


Equation (2), which is normally known as the biharmonic equation, may be 
written as 


The stress components, expressed in terms of a stress function, have the form 


Hence if, by definition, 


then the solution of equation (8) reduces to the solution of the following equations : — 


Oy ay 
axe 
¥ (11) 
and dy? = 
The finite difference form of these equations is given by 
(12) 
and | 


3. The Electrical Network 
3.1. ANALOGY 


The complete electrical resistance network is shown diagrammatically in Fig. 1. 
The low potential network consists of a square grid of side a, of equal resistances of 
value r. A network of this type alone would be suitable for solving the Laplace 
equation’. A similar network, but with resistances of value p, forms the medium 
potential grid. Corresponding points on each of these two grids are connected by 
large resistances of value R. In addition the high potential grid, which contains 
no resistances, has its nodal points connected to corresponding points on the 
medium potential grid by large resistances of value R’. The type of network formed 
by the R and r resistances alone, or alternatively the R’ and p resistances, would 
form a suitable analogy for the Poisson equation’. 
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Fig. 1. 
The electrical network, 


The complete network also makes provision for applying any desired electrical 
potential to the boundary points on both the low and medium potential nets, together 
with facilities for applying potentials to the nodes of the high potential grid. 

To show that this network does provide an analogy to problems of elasticity, it 
is only necessary to apply Kirchhoff’s law to two points in the network. 

Consider the total current flow into the node marked with potential V,, in the 
medium potential grid; then 
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from which 


+ 4 _ =0. . (14) 


It can now be seen that ince the three terms in brackets are small and can 
be neglected, then 


pE, 

and it will be realised that this equation is analogous to the first of equations (7) and, 
if E,=0, it is also analogous to the first of equations (12). 


VitVi+tV3+V, 


=0 ‘ (15) 


10 
6 2 5 
, 3 2 
| 
7 8 
le 
Fig. 2. 


Network pattern. 


Further, consider the total current flow into the node marked with potential v, 
in the low potential grid; then 


Here again, if the last term is negligible compared with the others, 


which provides an equation analogous to the second of equations (7) and the 
second of equations (12). 
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The precise physical functions which were taken to correspond to the electrical 
potentials are given in Table I together with the exact relationship between the 
physical and electrical parameters. 


If equation (3) is considered as it stands, instead of splitting it up into a pair of 
interrelated differential equations, then the finite difference form is 


4 8 12 qoa* 
=5 n=9 


n=1 n 


where the suffixes 1 to 12 refer to the value of w at points on a square grid of 
side a as shown in Fig. 2. 


The corresponding expression for v is given by 


(20) 


where the error is given immediately by the bracketed terms. 


3.2. CHOICE OF RESISTANCE VALUES AND ACCURACY OF ANALOGY 


The terms that are neglected in the electrical equations, numbers (15) and (18), 
can be made small provided that the resistance values are suitably chosen. Regard- 
ing equation (14), the fractional error can be seen to be of the order 


This can be made small provided that R and R’ are much greater than p. Since this 
condition also ensures that v is much smaller than V, then the second term will be 
particularly small, so that the fractional error will be of the order 


As regards equation (17), the fractional error is of the order 


r 


so that, provided that R is much greater than r, the error is truly negligible. 


For the present network the following values of the resistances were taken, 


r= 100 ohms, R=200,000 ohms 


p=100 ohms, R’=100,000 ohms and 
50,000 ohms 


These values ensured that the errors were negligible so that the analogies were 
truly justified. R’ was initially taken as 200,000 ohms but it was necessary to 
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TABLE I 
RELATIONSHIP BETWEEN PHYSICAL AND ELECTRICAL FUNCTIONS 
Flexure of Plate Extension of Plate 
Voltage 
Network Symbol || Order of — Order of Corre- 
Voltage | SPO"@!N8) Relationship || Voltage | 5P onding Relationship 
Plate Plate 
used used 
Function Function 
High 3R R’ 
Grid ial 
Medium MaR 2R 
Potential V M V ~ 0 
Grid 
Low 0->0-005 0->0-015 
Potential v w vu — volt 
Grid 


reduce this value, since the maximum potential available for the high potential 
net was 230 volts (E), and if R and R’ were too large, then v was too small to be 
accurately measured. 


It should be noted that the final effect of the resistances chosen for the network 
is to make the three potentials E, V and v to be all of different orders of magnitude. 


Table I gives the orders of the electrical potentials which were used in the 
experiments. 


3.3. CONSTRUCTION OF ELECTRICAL NETWORK 


That part of the network comprising the R, p and r resistances was built into a 
self-contained unit, a rear view of which is shown in Fig. 3. The mesh separation 
is one eighth so that there are 9 x9 nodal points. The terminals corresponding to 
the nodal points in the low and medium potential nets can be seen on the front of 
the unit in Fig. 4. This unit also contained some rheostats wired to the boundary 
points of the medium potential net to assist with the setting up of the boundary 
conditions. 


The R’ resistances comprised a separate unit which could be attached to the 
main unit when required. Potential dividers were also available for applying any 
required potential to the boundary points. All potential measurements were made 
with the usual potentiometer and galvanometer arrangement. The complete 
apparatus is shown in Fig. 4. High stability carbon resistors were used throughout, 
these being specially selected to give a tolerance of +4 per cent. about the 
mean values 
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Fig. 3. 
The electrical analogue. Rear view of main unit comprising the two principal 
nets and the interconnecting resistances. 


3.4. OPERATION OF ELECTRICAL NETWORK 
The operation of the network is simple and the procedure is as follows. 


If the problem is of flexure, then the potentials E must first be applied to the 
high potential grid by means of a potential divider; these potentials will be deter- 
mined from the intensity of loading by reference to Table I. If the problem is one 
of extension, the high potential grid is disconnected together with the R’ resistances. 


As regards the low and medium potential grids, it is necessary to set the 
potentials v and V, at the boundary points, to the appropriate boundary conditions. 
This involves a difficulty since it is not common for the boundary conditions of the 
problem to be in this convenient form. For example, instead of knowing v and V 
at the boundary, it is more likely that v and dv/0n, the normal derivative of v, will 
be specified. In this case, the procedure is to set the appropriate values of v on the 
boundary of the low potential net and then to adjust the values of V on the 
boundary of the medium potential net, until 0v/0n assumes the desired values round 
the low potential boundary. This involves a type of electrical iteration process, but 
with a little care and experience the desired conditions can rapidly be achieved. 


It is of interest to note that a similar electrical analogy has been proposed by 
Boscher and Malavard' and Liebmann™, but in both cases the method for setting 
up the boundary conditions was different. Liebmann fed currents into the low 
potential boundary points to set up the required slopes and then adjusted the 
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Fig. 4. 


The electrical analogue, showing all equipment (square plate with 
uniform normal pressure and built-in edges). 


medium potential boundary values to produce the required low potential boundary 
values; this method has also been tried by the present authors and is described 
later. Boscher and Malavard worked with zero boundary potential on the low 
potential grid, this being achieved by first producing a Laplace solution to an 
auxiliary problem. 


Although all problems investigated had two axes of symmetry, the board was 
not adapted for axes of symmetry as it was felt desirable to verify that identical 
results were obtained in each quadrant. 


4. The Experimental Results 
4.1. EXTENSION OF A PLATE BY A PARABOLIC DISTRIBUTION OF STRESS 


The particular problem of the extension of a flat plate, which was chosen to 
be investigated by means of the network, was that of a square flat plate subjected to 
a parabolic tensile stress distribution on each of two opposite edges. This problem, 
to which a theoretical solution exists, has been investigated analytically by strain 
energy by Timoshenko and Goodier“’, and by relaxation by Longman"?. 


If the square plate is of side d, with axes x and y from an origin at the centre 
of the plate, so that the loaded edges are at x= +d/2 and the unloaded edges are 
at y= +d/2, then the value of the stress function at the boundary is 


(2 
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Therefore, at the boundaries, the direct stresses are given by 


(c,) =0 


v= +4/2 


= (Tey) =0 


v=+d/2 z=4d/2 


The gradient of ¢ normal to the boundary is 


Ox/z=44/2 


(26) 


and x 


(3 v= 44/2 


Thus for the purpose of the tests, where v represents ¢, the constant K was taken 
to be equal to 0-12 volts, so that the maximum value of v on the boundary became 
0:01250 volts, and the slope normal to the boundary was zero on two edges and had 
the value 0:04/d=0-00500/a on the other two edges. 


To apply these conditions to the network, the high potential grid and the R’ 
resistances were removed, since they are not required in an extensional problem. 
The values of v on the boundary were then obtained from equation (24) together 
with Table I and applied, by means of the potential divider, to the low potential net. 
Finally, the values of V on the boundary of the medium potential net were adjusted 
until the corresponding values of dv/0n, on the low potential net, were as required. 


The setting of the boundary conditions was actually carried out in three different 
ways detailed as follows, the first two being similar to that already described and 
the third method similar to that proposed by Liebmann. 


(i) The first approximation to the gradient is given by 
Ov\ Up—V% 


where v, is the value of v on the boundary and vz is the value of v at the next node, 
a distance a from the boundary. This means that if v, and (0v/0n), are known, then 
vz is also known. Thus using this method v, can be set up on the low potential net 
boundary and the boundary potential on the medium potential net can be adjusted 
until vs assumes its prescribed value. 


The results of using this method are given in Fig. 5. 
(ii) A better approximation to the gradient 0v/0n is given by 


Ov\ 


24 The Aeronautical Quarterly 


) 
L 

2 
e 

|: 
| 
e 
tw 


ELECTRICAL ANALOGUE 


279 EXP |RES. 
376] THEORY 378] THEORY 


KEY KEY 


+ Vx 10° VoLT +Vx 10° Voir 


Fig. 5. Fig. 6. 
Flate plate with parabolic stress distribution on Flat plate with parabolic stress distribution on 
two opposite edges. (Gradient based on two two opposite edges. (Gradient based on three 
ordinates.) ordinates.) 


where vc; is the value of the node at a distance 2a from the boundary. Using this 
method, v, is set up as before, and then the boundary values on the medium 
potential net are adjusted until the relationship (28) is satisfied. This is rendered 
simple by means of an electrical circuit devised for the purpose and shown in Fig. 7. 
The slope, based on three ordinates, is given immediately by the potentiometer. 


The results of using this method are given in Fig. 6. 


(iii) A third possible method of applying the boundary conditions is to set the 
gradient of the potential at the boundary of the low potential net, before setting the 
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a | a J 
H/s H/, H af 
= ak (4-94 -4) 
oi 
AC (13) 1g 2a 1s 
H! 
Fig. 7. 
1 | Circuit for gradient based on three ordinates. 


boundary potential itself. This is achieved by feeding currents into the boundary 
points of the low potential net until they produce the desired gradient of potential. 
The boundary potentials on the low potential net are then obtained by adjusting the 
boundary potentials on the medium potential net. 


4.2. DISCUSSION OF RESULTS OF EXTENSION EXPERIMENTS 


The results of the experiments with the flat plate in extension are shown in 
Figs. 5 and 6, together with the corresponding theoretical values. Since the plate had 
two axes of symmetry the results corresponding to only one quarter of the plate 
were recorded, the other quadrants yielding almost identical values. The potentials 
on both nets were recorded, the actual voltages being given; the theoretical elasticity 
parameters being converted into corresponding electrical voltages for comparison. 
Also shown in the figures are the remainders or residuals obtained when the 
relaxation formula is applied to the experimental values. 


The expressions for the relaxation residuals are 


4 
n=l 
for the medium potential net, and 
ael 


for the low potential net. 
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In all three experiments, the relaxation residual is small, indicating that the 
analogy is accurate subject to the coarseness of the mesh used. In Fig. 5, where 
the gradient is based on two ordinates, the comparison between experiment and 
theory is not good, but in Fig. 6 where the gradient is based on three ordinates, 
the agreement is quite good in view of the small number of points used. In the 
case where the boundary conditions were applied in a different manner, similar to 
that used by Liebmann, the results were reasonable, but considerable practical 
difficulty was encountered since there was a large degree of cross-coupling on the 
boundary conditions, such that when the boundary potential was adjusted the slope 
was also affected and vice versa. Also, it was noted that the potentials were less 
steady, due to the boundary circuit arrangement, and this caused the finite difference 
formulae residuals to be higher. This method of setting the boundary conditions 
was not proceeded with. 


4.3. FLEXURE OF A PLATE BY A UNIFORM DISTRIBUTED LOAD 

The particular problem of the flexure of a plate which was investigated was that 
of a square flat plate subjected to a uniform normal loading, the edges of the plate 
being either simply-supported or built-in. Again, the theoretical solution of this 
problem is known and is given by Timoshenko", 

In these experiments the high potential grid was required, together with the R’ 
resistances. A uniform distributed load required all the nodes in this grid to be 
connected together and a large potential applied to them, since E is to be constant. 

Dealing first with the case of simply-supported edges, the relevant boundary 
conditions are 


M,=M,=0 (31) 
and w=0 
so that as regards the electrical network, the conditions become 
(32) 


These conditions are particularly simple to apply as they imply that the 
boundaries of both the medium and low potential nets are shorted and earthed. 


This experiment was carried out twice, once with R’= 100,000 ohms and once 
with R’=50,000 ohms, so as to increase the values of the potentials v. The results 
of the latter experiment are given in Fig. 8. 


For built-in edges, the boundary conditions are 
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Fig. 8. Fig. 9. 


Square plate with uniform normal pressure. 


Simply-supported edges. 


Square plate with uniform normal pressure. 


Built-in edges. (Gradient based on three ordinates.) 


so that for the electrical network the conditions become 


and 


n 


(34) 


These conditions were applied in the same way as were the boundary conditions 
for the case of a plate in extension; Fig. 9 gives the results of the experiment when 
the boundary gradient is based on three ordinates, with R’=50,000 ohms. 
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4.4. DISCUSSION OF RESULTS OF FLEXURAL EXPERIMENTS 


The results of the experiments with a flat plate in flexure are given in Figs. 8 
and 9. Here, owing to symmetry, only one eighth part of the nets were recorded, 
the other sectors yielding almost identical results. The finite difference formulae 
residuals, together with the theoretical results, are also given, all in terms of 
electrical voltages. 


For simply-supported edges the results are extremely good, but in this case it 
was not necessary to perform any electrical iteration, since it was possible to apply 


directly the boundary conditions. A 


The results for built-in edges were not so good, even when the gradient was 
based on three ordinates. However, this is thought to be mainly due to the coarse 
net used, since only nine ordinates are available to define a curve with two points 
of contraflexure. 


It was also noted that the high potential E, which was a mains voltage, tended 
to drift and so make the recording a little tedious; this voltage should, for preference, 
be stabilised. 


It should be noted that if the boundary conditions specified involve a differential 
higher than the first, then a procedure similar to that employed for setting a gradient 
based on three ordinates could be used. This condition would apply if the plate 
had a free edge, but this experiment was not tried as it was felt that the network 
available was too small to be likely to yield accurate results. 


5. Conclusions 

The work detailed in this report shows that the use of a pure resistance net- 
work, to solve some of the problems of extension and flexure of flat plates, is quite 
practicable. 


The results obtained are quite satisfactory, bearing in mind the size of the 
network used. The fact that the finite difference formulae residuals are small shows 
that the solutions obtained with the network will be practically identical with those 
obtained by a relaxation solution on a grid of the same size. 


The network is simple to construct and, since it only consists of resistances, it 
is very reliable. The operation of the network is a simple matter. 


If a large network were to be built it would be desirable to have voltage 
stabilisation for the high potential net. Also, if the number of boundary points were 
increased it might be necessary to have a display of the boundary potentials, possibly 
by means of a battery of millivoltmeters. 
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On an Approximate Method of Calculating 
Pressures on Non-lifting Head Shapes at 


Supersonic Speeds 


H. K. ZIENKIEWICZ 
(Fluid Motion Laboratory, University of Manchester) 


Summary: Slender-body theory is used to derive the ogive of curvature 
approximation for very slender, pointed, convex head shapes at supersonic 
speeds. Results of application of this approximation, together with the 
A-method for circular arc ogives, to a variety of non-slender head shapes 
show very good agreement with the method of characteristics, van Dyke’s 
second-order theory and experiment. Good agreement with the method of 
characteristics and with experiment is obtained even in cases when the 
stagnation pressure losses across the nose shock wave are not: negligible. 


1. Introduction 


In Ref. 1 a number of approximate methods of calculating pressures on non- 
lifting, pointed, convex head shapes at supersonic speeds have been put forward. 
The present paper summarises some recent developments of the most promising of 
these methods, the so-called “ ogive of curvature method.” In particular, the basic 
assumption underlying this method is justified by showing that it leads to the exact 
result in the limit of vanishingly small supersonic similarity parameter ft, the 
error in the pressure coefficient at any point on the head relative to the pressure 
coefficient at the nose vanishing as (log 8t)-'. This result, obtained by using the 
slender-body theory, is rather more satisfactory than the original justification of 
the assumption, which was based on purely empirical grounds". 


The ogive of curvature approximation is then used together with the A-method 
for circular arc ogives®’ to calculate pressures on a number of non-slender head 
shapes at Mach numbers such that the supersonic similarity parameter ft varies 
between about 0-4 and 1:0. Comparison with the corresponding pressure distribu- 
tions calculated by the method of characteristics (taking account of flow rotation) 
and by van Dyke’s second-order theory and with some experimental pressure 
distributions shows very good agreement and leads to the conclusion that the 
ogive of curvature method effectively bridges the gap between the linearised and 
the conical-shock-expansion theories’. 


NOTATION 
A,B,....,H designation of head profiles (see Table I) 
C, pressure coefficient, =(p—p,)/q 
Cy.  ogive of curvature pressure coefficient 
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K = 28tan(x/2), supersonic similarity parameter for circular arc 
ogives 


K, = 2M tan(x/2), hypersonic similarity parameter for circular arc 
ogives 


L _iength of ogive of curvature 
M free stream Mach number 
p.p, local and free stream static pressures, respectively 
q dynamic head, =4yp,M? 
r(x) equation of head profile 
s,S cross-sectional areas of head and ogive of curvature, respectively 


t,T thickness ratios (maximum diameter/length) of head and ogive 
of curvature, respectively 


x,randX,R__ cylindrical co-ordinates (see Fig. 1) 
y,Y  x-and X-like variables of integration 


B = /(M?-1) 
y ratio of specific heats ‘ 
6 = tan-'?’ (x), profile slope t 


X factor of proportionality defined by equation (9) 
x  ogive of curvature nose semi-angle 
Dashes denote differentiation with respect to x or X. 


Suffixes 
t two-dimensional simple-wave flow 


n conditions at the nose 


©) 


2. The Ogive of Curvature Approximation in the Slender- 
Body Theory 


Referring to Fig. 1 consider a slender, pointed, convex head profile at zero 
incidence in a supersonic air stream of Mach number M. The length of the body is 
taken as unity and its maximum diameter is t. The equation of the head profile is 
r=r (x), such that 7’ and r” are O(f) and it is assumed that 7’ and r” are everywhere v 


continuous but r’ (1) is not necessarily zero. Then, under the usual assumptions of a. 
the slender-body theory, the pressure coefficient at a point P(x,r) on the body te) 
surface is given by“? e 
a 
C= Glog x+s” (x) log + | og x-yats” onl} - 
0 

. ‘ (1) B 
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Fig. 1. 
Co-ordinates and notation. 


where s=ar* is the cross-sectional area of the head shape and 8=1/(M?-1). In 
equation (1) the dominant term is the second, being O (t* log t), while the remaining 
terms are O(t?). Hence, we can rewrite equation (1) as 


— 


Consider next the flow about the ogive of curvature of the head profile at P 
(Fig. 1). This ogive is obtained by revolving about the axis of symmetry an arc 
of the circle of curvature of the profile at P. The length of the ogive up to the 
point of zero slope is L and its maximum diameter is TL. Approximating the ogive 
of curvature by the corresponding parabolic arc body (differences in ordinate, slope 
and curvature between the two profiles are O(T*) and therefore are negligible in 
the present application), it can be shown that 


T=} 


where the orders of magnitude follow from the assumption that r, r’, and r”’ are 
all O(#). Thus, the ogive of curvature is also a slender body so that, with the origin 
of the X, R co-ordinates at the nose of the ogive and approximating as in 
equation (2), the pressure coefficient at the point P(X, R) on the ogive of curvature 
at P is 


Ss” 


Xi0om. . . .. @ 


BR 


But, by the definition of the ogive of curvature, R(X)=r(x), R’(X)=r (x) and 
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R”’ (X)=r’ (x), so that S’(X)=s’(x). Hence, by virtue of equations (3) it 
immediately follows from equations (2) and (4) that 


C,=Cy. + O (t?) ‘ ‘ (5) 
provided that r’ and r’ do not vanish simultaneously, for then X =00. 


Thus, to order f? log t the pressure coefficient at a point P on a slender, pointed, 
convex head shape is the same as that at P on the ogive of curvature at P at the 
same free stream Mach number. The error relative to the pressure coefficient at 
the nose of the body, which is proportional to flog Bt, vanishes as (log Bt)-* and 
the ogive of curvature approximation becomes exact in the limit of 8t—> 0. 


At the nose, the actual profile and its circle of curvature have the same slope 
and curvature, and near the nose they coincide to first order in x. Hence, obviously, 
the ogive of curvature approximation gives the correct pressure coefficient at the 
nose. Further, since for a given Mach number the pressure gradient at the nose 
depends only on r (0) and r” (0), provided that the flow is everywhere supersonic 
(see, e.g., Ref. 5 for a full discussion), the ogive of curvature approximation gives 
also the correct value of this gradient. Hence equation (5) can be written more 
precisely as 


+O f@ri . ‘ ‘ (6) 
where f(x)=O(1) and depends on the profile shape. 


In fact, the ogive of curvature approximation in the slender-body theory gives 
better accuracy than might be expected from the estimate of the order of error, 
because some of the O(r*) terms cancel out (e.g., the r’? (x) term in equation (1), 
although of O(t?), is given correctly by the ogive of curvature approximation, 
because r’=R’). It appears difficult to obtain an estimate of the error for an 
arbitrary profile better than that in equation (6), due to the awkward form of 
equation (1) and of the relation between x and X (equations (3)). However, the 
following example may serve as an indication of the accuracy of the ogive of 
curvature approximation in the slender-body theory. 


Take the equation of the head profile to be given by r=4tx(3—2x)!. Then, 
from equation (1) 


while, applying the ogive of curvature approximation, it is found that 
C,.=C,+ AC, 
where C, is as in equation (7), and 


AC, 3 (1 — x) [3 (3 


3(1—x)+[3 (3—2x+ 
2 (3 — 2x) ‘ 


2 (3 — 2x) 


— (1 —2x) log 


(8) 
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Bt=0°1 


A 


— SLENDER-BODY THEORY 
EQUATION(8) 


0-2 0-4 0-6 0-8 1:0 
x 
Fig. 2. 
Comparison of pressure distributions as given by the slender-body theory 
and the ogive of curvature approximation, 


C, and C,, for Bt=0-1 are plotted against x in Fig. 2. It can be seen that over the 
forward 80 per cent. of the head length AC, is less than 5 per cent. of C, at the 
nose; the relative differences between C, and C,, can, of course, be made arbitrarily 
small by decreasing ft. 


3. Application of the Ogive of Curvature Method to 
Non-Slender Head Shapes 


The ogive of curvature approximation is of little practical value in the slender- 
body theory, except perhaps for head shapes with s” such that the integral in 
equation (1) cannot be evaluated in a closed form. What merit the method 
possesses lies largely in the fact that it can be applied to non-slender head shapes 
for which no other equally rapid and accurate method exists at present. The 
application of this method to such head shapes is made possible by the availability 
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0-4 
0 
0 0:4 0-8 1-2 1-6 2:0 
K= 2B tan (x/2) 
Fig. 3. 


Circular arc ogives. Variation of A with supersonic similarity 
parameter and free stream Mach number. 


of accurate pressure distributions on non-slender circular arc ogives. One could 
use directly the characteristics pressure distributions, generalised on the basis of 
the hypersonic similarity law, which are given in Refs. 6 and 7. However, since 
the similarity generalisation involves a certain loss of accuracy, particularly at 
Mach numbers between 1:5 and 2, it is more accurate to use the so-called A-method 
put forward by the author’. For convenience of readers unfamiliar with this 
method a brief summary is given in what follows. 


On analysing the characteristics pressure distributions of Ref. 6, the author 
discovered’? a remarkably simple relationship between the axi-symmetric flow 
over a circular arc ogive and the two-dimensional simple-wave flow over the 
same profile. Let p, and M, be the static pressure and the local Mach number, 
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0°30 
—— |CHARACTERISTICS 
|\—METHOD 


0-25 


0-20 \ 


Cp 
0:10 
N 
0:05 
\ 
= 
—0:055 O-4 0-6 0-8 
Fig. 4. 


Circular arc ogives. Typical comparison of A-method with 
method of characteristics. y=18-9°, M=3. 
respectively, at the nose of a circular arc ogive, p be the static pressure at a point P 
on the ogive in axi-symmetric flow and p, be the static pressure that would obtain at 
the point P if the flow expanded two-dimensionally over the same profile from the 
given conditions at the nose. Then, it turns out that 


Pu— Pr) 
or Co =A Cy) 


where the factor of proportionality A is constant over the profile and a function only 
of the free stream Mach number M and of x, the nose semi-angle of the ogive. 


(9) 


In Ref. 2, A was calculated from the irrotational characteristics results and a 
chart was given showing the variation of A with x and M for 15=<=MX3:°5, 
0-5 = 2M tan (x/2)=Kn<1-0. This chart has since been revised and extended to 
K,=2 by using the characteristics results of Ref. 7, which included the effects of 
flow rotation, and to values of K;, <0-°5 by using a number of unpublished charac- 
teristics and second-order theory solutions. The revised and extended A chart is 
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shown in Fig. 3; the supersonic similarity parameter K, rather than the hypersonic 
parameter K,, is used as the abscissa to reduce the dependence of A on M. The 
value of A for K=0-1 was obtained using the slender-body theory. Above this 
value of K the slender-body theory is generally considered to be unreliable, at least 
as regards pressure distributions. For K much below 0-1 this theory indicates that 
A ceases to be approximately constant along the profile. Thus K=0-1 appears 
to be the lower limit of the applicability of the A-method. In addition, the method 
is at present limited to M> 1-5 owing to the lack of data required for evaluating A 
at M< 1-5. However, it is believed that equation (9) would still apply at Mach 
numbers below 1-5 provided that x is small enough for the flow over the ogive to be 
everywhere supersonic and large enough for K>0-1. At the other end of the range 
of K for which A has been determined, it is of interest to note that as K increases 
above about 1-5, A very closely approaches unity, indicating that the flow is then 
practically two-dimensional; this is in agreement with Eggers’s conical-shock- 
expansion theory for high Mach number flows". 


As an indication of the accuracy of the A-method, a typical comparison with the 
method of characteristics is shown in Fig. 4. 


The procedure for applying the A-method, together with the ogive of curvature 
approximation, to calculation of pressures on non-slender head shapes is very simple. 
Once the distributions of the profile slope 


6=tan-'?r’ (x) 
and of the ogive of curvature nose semi-angle 
x=cos~? [cos 6(1+rr’” cos? 


are known along the profile, the corresponding values of A for a given Mach number 
are found from Fig. 3. Then, using the tables of flow about cones‘? to find 
Cyn (x, M) and a standard Prandtl-Meyer expansion table to find C,, by the usual 
methods, C,, the pressure coefficient at any point on the profile follows from 
equation (9) and the ogive of curvature approximation. Computing details, suitable 
pro formae and worked examples are given in Ref. 1. 


When applying the ogive of curvature method to head shapes such that the loss 
of stagnation pressure across the nose shock is not negligible (greater than, say, 
5 per cent. of the free stream stagnation pressure) one of two alternative procedures 
can be followed. The first is to take the stagnation pressure on the local ogive of 
curvature as equal to that at the nose of the actual body, which requires different 
free stream stagnation pressures for the actual body and the ogive of curvature; the 
second is to take the free stream stagnation pressures as equal, which violates 
the condition of constant stagnation pressure along the body surface. At first sight 
the first alternative appears more plausible and was in fact adopted in Ref. 1, since 
for all the head shapes considered at the time stagnation pressure losses were 
negligible and there was no evidence as to which procedure would lead to better 
accuracy. However, since then investigation of pressure distributions on head 
shapes with large nose angles (and at Mach numbers such that losses of stagnation 
pressure across the nose shock were up to 30 per cent. of the free stream stagnation 
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TABLE I 
PROFILES OF HEAD SHAPES 
Fineness 
D 
esignation Equation Ratio (1/1) 
r= 3 
B 3 
0-03968 (7x + sin 7x) 4+ 
+ 

D @—4sin2y 1-22) 3 
E [x (2—x)}*/4 3 
F 3 
G 0:0655 x + 0:02782 logio (1 + 42°37 x) 45 
H 0:0636 x + 0:03286 logio (1 + 26°74 x) 45 


pressure) have revealed that the first procedure consistently underestimated pressures 
on such heads and that the second procedure yielded much better accuracy. Con- 
sequently, the second procedure has been adopted as standard and in fact is 
implicit in the ogive of curvature approximation as stated in the present paper, 
namely, that “the pressure coefficient at a point P on a convex head shape is the 
same as that at P on the ogive of curvature at P at the same free stream Mach 
number”. This modification of the original method requires an obvious modification 
of the computing procedure given in Ref. 1. If velocity or Mach number distri- 
butions along the head are required, these should be obtained from the pressure 
distributions using the correct value of the (constant) stagnation pressure along the 
body surface and not the fictitious values on the local ogives of curvature used in 
the computation of pressure coefficients. 


A difficulty arises in applying the ogive of curvature method to head shapes 
such that curvature of the profile becomes infinite or zero at the end of the head*. 
In the former case, —r” (1)=00, as x—>1 a point is reached where x is so large 
that it falls outside the range of the A-y chart (Fig. 3) and even if this were not so, 
curvature would eventually increase sufficiently for the whole concept of the ogive 
of curvature to become meaningless. When r” (1)=r’ (1)=0, x —>0 as the end of 
the head is approached and a point is reached where K <0:1, which is outside the 
limit of applicability of the A-method. In the Appendix, a procedure for dealing 
with such cases is suggested which, though perhaps somewhat artificial, appears to 
give acceptable results. In any case, the region of a head affected by —r”’ (1)=0 
or 00 is usually small and it would not matter very much how one extrapolated there 


*Cases when this occurs elsewhere on the head need not be considered here. For, profiles with 
infinite curvature at 0<x<1 are unlikely to be of practical interest, the possibility of 
r’(0)= — oo is excluded by the assumption of a sharp nose, while the case of zero curvature 
presents no difficulty unless the slope is also zero (cf. equations (5) and (3)), which on a convex 
profile with continuous slope can only occur at x=1. 
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0-40 
| 
o-35|-4 
—--]| OGIVE OF CURVATURE METHOD 
FROM REF. 3 FOR 
___ | CHARACTERISTICS ] 
0-30 
0-25 
0-20 
C, \ 
0-15 
PROFILE E 
0-10} 
\\ 
N 
0-05 — ee 
= S J 
PROFILE H 
“0°05, 0-2 0-4 0-6 0-8 
Fig. 5. 


Comparison of ogive of curvature method with method of characteristics. 
Profile E at M=3-0 and profile H at M=2°3. 


the pressure distribution, particularly from the point of view of wave drag calcula- 


tions (the profile slope near the end of the head being necessarily small, especially 
when r” (1)=r’ (1)=0). 


4. Comparison with More Exact Theories and Experiment 

Since it was first put forward, the ogive of curvature method has been extensively 
tested by comparison with results of other theories and with experimental results 
on several head shapes. Since in the ogive of curvature method the pressure 
distribution is closely connected with the profile slope and curvature, profiles with 
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OGIVE OF CURVATURE METHOD 


0-28-74 

|| ——]SECOND-ORDER THEORY [FROM REF.1!] 
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24 0:24 
— SECOND-ORDER THEORY 
—- OGIVE OF CURVATURE METHOD \ (b) 
0-20 0-20 
\ 
\ 
q 
\ 
() \ 
0-16 
WN 0-16 
\ \ 
0-12 \ A \ 


0-12 
PROFILE A, PROFILE F, M=I-5 
Cp V4 Mz2 


Ww PROFILE F, M=3 


\ [AC prorive 8 
N 
PROFILE 
3 M=2 
0 NS i\ 


| \ 
-0-04 


-0:08 


0-2 0-6 08 ro ~0°08 
O-2 0-4 0-8 1-0 
x 
Fig. 6. 
Comparison of ogive of curvature method with second-order theory. 
(a) Profiles A, B and C. (b) Profiles D and F. 


a wide variety of slope and curvature distributions were chosen for comparison. For 
obvious reasons it is possible to include here only a representative selection of the 
results. Profile geometry of the head shapes considered is summarised in Table I. 


4.1. METHOD OF CHARACTERISTICS 

Figure 5 shows pressure distributions obtained by the ogive of curvature 
method on profiles E and H compared with the corresponding pressure distributions 
calculated by the method of characteristics, including the effects of flow rotation. 
The profile E has a blunt nose and, to make it possible to apply the method of 
characteristics, a cone tangent to the profile at x=0-05 was used in Ref. 3 to replace 
the blunt nose. When applying the ogive of curvature method it is, of course, more 
convenient to replace the blunt nose by a circular arc ogive rather than a cone; this 
was not done in the case of the profile E, since the pressure distribution for x < 0-05 
is of no interest for the present purpose. 


Particularly gratifying is the good agreement obtained for the profile H. This 
profile has a large nose semi-angle (24°) followed by a very rapid change of curvature 
near the nose, while the curvature and slope near the end of the head are so small 
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0:28 
0-24 
— |OGIVE OF CURVATURE METHOD 
[REF 13] 

0-20} 

O-16 

Cp 

0-12 

0-08 \ 

0-04 

0 
© 0-2 0-4 1-0 
x 

Fig. 7. 


Comparison of ogive of curvature method with experiment. Profile G at M=2-45. 


that the limit of applicability of the A-method is approached there. Hence, com- 
parison of pressure distributions on this profile is a severe test of the ogive of 
curvature method. Comparison of the computing labour required by the two 
methods is also of interest. Calculations by the ogive of curvature method took 
approximately 1 hour. In the characteristics solution* a very small mesh size had 
to be used near the nose, which resulted in some 280 mesh points and a computing 
time of over 300 hours. 


4.2. SECOND-ORDER THEORY 


By now, van Dyke’s second-order theory"® has been widely applied by several 
investigators and, in the case of head shapes, has been shown to give very good 
accuracy for 8t<1. Although much more rapid than the method of characteristics, 


*Carried out at the Mathematics Division, National Physical Laboratory, under the direction 
of A. R. Curtis. 
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it still requires a considerable computing effort, of the order of 10-20 hours for one 
pressure distribution. 


Good agreement between the ogive of curvature method and the second-order 
theory has been obtained in all the several cases investigated, some of which are 
shown in Fig. 6(a) (profiles A, B and C) and Fig. 6(5) (profiles D and F). Of special 
interest are regions near the end of the heads A and D, which terminate with infinite 
curvature of the profile, and of the head C, which terminates with zero profile slope 
and curvature. In these regions the procedure referred to in the last paragraph of 


Section 3 and described in the Appendix was used with apparent success. 


4.3. EXPERIMENT 


The accuracy of the ogive of curvature method has also been explored by 
comparison with wind tunnel experiments on a number of head shapes. The 
results, published in Refs. 12 and 13, show good agreement. A representative com- 
parison (profile G at M=2:-45) is shown in Fig. 7. The profile G is similar to the 
profile H (cf. Table I), but has an even larger nose semi-angle (30°) and the free 
stream Mach number is such that the stagnation pressure loss across the nose shock 
is about 14 per cent. of the free stream stagnation pressure. As can be seen from 
Fig. 7, the theory predicts correctly the enormous expansion from the pressure at 
the nose (C,,=0-600) and over most of the head agreement with experiment is 
within the experimental scatter*. Similar agreement was obtained for other head 
shapes and Mach numbers"®?. 


5. Conclusions 

Slender-body theory shows that the ogive of curvature approximation becomes 
exact in the limiting case of head shapes with vanishingly small supersonic similarity 
parameter St. When this approximation is used together with the A-method for 
circular arc ogives, comparisons with the method of characteristics, van Dyke’s 
second-order theory and experiment show that its range of applicability is much wider 
than might be expected from the assumptions made in its derivation and that in fact it 
gives very good accuracy for non-slender head shapes with ft as high as 1-0-1°5. 
Since for Bt> 1-5 Eggers’s conical-shock-expansion theory becomes accurate®:®, 
the ogive of curvature method bridges the gap between Eggers’s theory and the 
linearised and slender-body theories. 


Interpreted as an extension of the familiar concept of two-dimensional simple- 
wave flow, the ogive of curvature method implies that the pressure coefficient on a 
convex axi-symmetric head shape depends only on the free stream Mach number, 
the local profile slope and the local ratio of the radius of curvature of the profile to 
the profile ordinate. 


In spite of its accuracy the method is very rapid, requiring only about a tenth 
or a hundredth of the computing effort necessary for the application of the second- 
order theory or the method of characteristics, respectively. 


*Four pressure measurements were obtained at each station as the model was rotated through 
four angular positions®. In Fig. 7 the average of these readings is indicated by a dot and 
their scatter by a vertical stroke. 
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There is nothing in the slender-body theory derivation of the ogive of curvature 


approximation to suggest that it should not apply to complete bodies of revolution 
as well as head shapes, provided that profile slope and curvature are everywhere 
continuous. However, owing to the lack of accurate pressure distributions beyond 
the point of maximum thickness on circular arc ogives, it is not possible to say at 
the moment whether this would also hold for non-slender bodies. 
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Appendix 


TREATMENT OF PROFILES WITH ZERO OR INFINITE CURVATURE 
AT THE END OF THE HEAD 


Referring to Fig. 1, let the point P now denote the last point on the head for which 
A can be determined from Fig. 3; at x=1 we have either r” = — 00, or r’=r’=0. 


Differentiating equation (9) with respect to 6 yields 
dp dn. 
«dé 


Using the well-known Prandtl-Meyer relation this can be written as 


M, and M, are the Mach numbers corresponding to p, and p, respectively, and 


(10) 


dp ypM;? 
d6/, /(M,?-1) 


is the local pressure gradient, with respect to flow direction, which would exist in the 
two-dimensional simple-wave flow with local static pressure p and Mach number M,. 


Thus, equation (10) gives the axi-symmetric pressure gradient at the point Pasa 
fraction, A’, of the corresponding two-dimensional pressure gradient. 


When r”(1)= — 00, the flow is locally two-dimensional at x=1 (this follows from 
the characteristics equations) and A’=1 there. Since A’ is known at P and can be 
similarly found for any point upstream of P, A’ can be plotted against @ and its average 
value, A’,, between @p and (@),., found, e.g. by graphical integration. Then the change 
of pressure between the point P and the end of the head is given approximately by A’, 
times the corresponding two-dimensional change of pressure. 


When r’(1)=r(1)=0, x—>O as x—>1 and for values of yx such that 
28 tan (x/2) < 0-1, A is no longer constant along local ogives of curvature (cf. Section 3). 
However, the average value of A tends to zero as y—>0. This follows simply from the 
fact that for small , the axi-symmetric change of pressure between the nose and the end 
of a circular arc ogive is proportional to x? log By, while the corresponding two- 
dimensional change of pressure is proportional to x; hence the ratio of the two changes 
of pressure, which is a measure of the average A, varies as x log By. Therefore, we 
can proceed as in the case of r”(1)= — 00, but now we take A’=A=0 at x=1. 
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Some Properties of Nodes in Vibrating 


Systems 


R. E. D. BISHOP, M.A., M.Sc.(Eng.), Ph.D. and D. C. JOHNSON, M.A. 
(University Engineering Laboratory, Cambridge) 


Summary: Special forms of the receptance functions of mechanical systems 
are shown to indicate the existence of nodes. These are associated with 
horizontal straight lines in receptance diagrams. The effects of this nodal 
condition are explained in terms of simple torsional and translational 
systems for the sake of simplicity; but the theory may be extended to any 
type of system and any set of generalised co-ordinates. A particular instance 
arises if principal co-ordinates are used. _, The paper is a sequel to an earlier 
paper®), 


1. Introduction 


Receptance* functions have been used by several writers®:*?. They have 
been shown to lead to great simplification of vibration analysis when complicated 
systems are concerned. The purpose of this paper is to explain the special, 
degenerate analytical forms that receptance functions may take and to show how 
their existence is manifested in receptance diagrams. The technical importance of 
this is explained in Section 3. 


It is possible to develop the theory, which is given here, in terms of any set of 
generalised co-ordinates; a case of special interest arises if principal co-ordinates 
are used. In fact, this has been done by the authors in a book which will appear 
in due course. This plan is not followed here, partly for the sake of simplicity and 
partly because the main features of the theory can be given with reference to 
simple torsional oscillations of conservative systems with finite freedom. For this 
purpose, results can conveniently be taken from Bishop’s paper” of which the 
present paper can be regarded as a sequel. 


NOTATION 
A,B,C.... vibrating systems 
xA,, constant defined in equation (6) 
F, amplitude of torque applied to s” disc 
f (w’), F (w’), etc. polynomial functions of w? defined as required 
polar moment of inertia 


*The word “receptance” is used instead of ‘“ mechanical admittance” in accordance with a 
suggestion which has been published; e.g. see Journal of the Royal Aeronautical Society, 
April 1954. 


Received June 1954. 
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k spring stiffness (see Figs. 4 and 5) 
M _ mass of particle (see Figs. 4 and 5) 
p, r™ principal co-ordinate 
t time 
x, displacement of s“ disc (also see Fig. 4) 
a,B,y,... receptances of A, B,C.... 
4 polynomial in w? defined by equation (3) 
A spring stiffness (see Figs. 4 and\5) 
excitation frequency 


2. The Series Form of Receptances 


If rigid discs are attached to a light elastic shaft so as to form a simple 
torsional system which is free from damping, that system may be thrown into 
oscillation by a harmonic torque F,sinwt which is applied to the s“ disc. If 
attention is confined to the steady harmonic oscillation having the applied frequency 
w, then the motion of the r™ disc is found to be proportional to F,sinwt and may 
therefore be expressed by the equation 


where 2,, is a constant called the “cross-receptance” between x, and x,. Its value 
is found to depend upon the constants of the system concerned and upon o. If the 
motion is measured at the same disc as that to which the torque is applied (the r™ 
disc, say) then 


where 2,, is now the “ direct-receptance” at x,. 


Table I of Ref. 2 contains the receptance functions of several simple torsional 
systems. It is seen that all are of the form 


where f (w) and A are polynomials in w. It is also seen that all the receptances of 
any given system have the same denominator 4. The frequency equation of the 
system is given by 

A=0 ‘ ‘ (4) 


and it will be supposed to have roots w?=w,’, ,’,....,’; there are as many of 
these as there are degrees of freedom—that is, there are n of them—and they are all 
real and positive. It follows that equation (3) may be written 


fw’) 
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Finally, it is found that, as any receptance of a system can be written in this form, it 
can be expanded into a series of partial fractions such that 


. 


where the numerators are all constants. The quantities ,, ,,....@, are the 
natural frequencies of the system in order of ascending magnitudes. 


These observations can, of course, be generalised. This can be done by means 
of Lagrangean theory (see Ref. 4) but, for the purposes of this paper, the results can 
conveniently be quoted. 


Suppose that a system of n discs mounted on a light elastic shaft is brought into 
steady harmonic oscillation by a harmonic torque F,sint applied at the s disc. 


The angular motions x,, x.,...Xn of the discs are given by 


x,=4,,F,sinwt 


X,= %,F, sin wt (7) 


Xn=%;F, sin ot 


By the series form of the receptances, therefore 


8 8 
n=[ A, + 2A +....+ | Fesinot 


x,= [= = F.sin ot 


If, now, » approaches the k™ natural frequency w,, the k column of fractions will 
account for more and more of the motion and, in order that the motion shall remain 
finite, the amplitude F, of the applied torque must be diminished. By this limiting 
process, a state of free oscillation with frequency , is approached such that 


That is to say the relative magnitudes of .A,., .A2s,...~Ans define the shape of the 
k™ principal mode. It may be mentioned that the same modal shape is reached no 
matter which disc is acted on by the harmonic applied torque. That is, the 
subscript s of equation (9) plays no part in defining the modal shape. 
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3. Receptance Diagrams 


It is convenient to discuss two types of graphical representation of direct recep- 
tances. The first is a plot of z,, against w and the second is a plot of 1/2,, against o. 


A typical curve of the first kind is shown in Fig. 1 in which © is plotted 
vertically. The full-line branches of this curve are found by adding the dotted 
curves in accordance with equation (6) when the equation is that of a direct 
receptance 2,,; each dotted curve represents one of the partial fractions and thus 
has the form of the resonance curve of a simple oscillator. Fig. 1 illustrates the fact 
that the constants ,A,,, .A;r,...nArr all have the same sign. The reason for this is 
that, were it not so, loops would appear in the diagram which would mean that one 
or more branches was multi-valued in <,, or w; this cannot be so for a reason which 
will be explained later. 


The second important type of receptance diagram is that of 1/2,, against w, and 
a typical curve is shown in Fig. 2. This curve is found by plotting horizontally the 
reciprocals of the abscissae of the curve of Fig. 1. The natural frequencies of the 
system are given by the n intersections of the branches of this curve with the w-axis, 
where ¢,, is infinite. On the other hand, the n—1 intermediate frequencies at which 
there are horizontal asymptotes are “ anti-resonance frequencies ”; at these values of 
©, %,=0, so that a harmonic torque applied to the r™ disc (with one of thse 
frequencies) will produce no motion of that disc. 


This second type of diagram has been used extensively (e.g. see Refs. 1, 4 and 2). 
It has also found application in certain numerical techniques. The reason for 
the importance of diagrams of 1/2,, is that they are helpful in the discussion of 
complex systems. Thus, if a system A can be considered as made up of two 
sub-systems B and C which are linked at a single co-ordinate x, (say), then 


1 
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where the receptances of A are denoted by «, those of B by £ and those of C by y. 
Therefore, the abscissae of the 1/8,, and 1/+y,, diagrams may be added to give the 
1/a,, diagram*. 


Of particular interest are those cases where B is some given system and C is a 
particle, an anchored spring or a simple oscillator. For instance, by contemplating 
the effects of making such modifications as these to any vibrating system B, it is 
easily found that the branches of Fig. 2 cannot be multi-valued in 1/2,, or o (6,, 
being identified with <,, in this case); this is therefore true also of the branches of 
Fig. 1. Again, if the free end of an anchored massless spring of stiffness k is 
attached to B at x,, y,-=1/k and the frequencies of the modified system are given by 
the equation 


1 
Bre 


The full lines in Fig. 2 now represent branches of the (1/8,,, o) curve and the curve 
for 1/7, is the vertical dotted line. It is seen that all the natural frequencies of the 
original system are raised by the presence of the spring. The addition of this spring 
will be referred to as the “addition of stiffness” at the co-ordinate concerned. 


The shapes of the receptance curves for a given system are determined by the 
constants of the system and may vary within the limits mentioned. Now it is 
often useful to know the shapes of these curves as may be seen from the following 
simple example. Suppose that the full-line curve of Fig. 2 represents 1/8,,, the 
reciprocal of the direct receptance at the r disc of a torsional system B, and it is 
required to estimate the effect upon the natural frequencies of B of increasing the 
polar moment of inertia of the r" disc, by an amount J. This increase of inertia can 
be regarded as a system C such that 


1 1 
Yrr 


According to equation (10), the natural frequencies of the modified system are 
given by 


1 
Bee 


That is, the intersections of the full-line curves and the dotted parabola in Fig. 2 
determine the new natural frequencies. Now, if a particular full-line branch is 
sensibly flat for positive values of 1/8,, and only curves upwards to the left of the 
axist, then the natural frequency appropriate to the intersection of that branch with 
the vertical axis is not altered materially by the attachment of the extra inertia (that 
is, by system C). 


*This feature is dealt with at length in Ref. 2. 
+This is the case with the upper branches of Fig. 11 of Ref. 2. 
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NODES IN VIBRATING SYSTEMS 


It is natural to enquire to what extreme shapes practical receptance curves can 
tend. The rest of this paper is devoted to a discussion of these special forms of the 
curves and to the implications of the underlying theory. 


The diagrams mentioned are concerned only with direct receptances. It is 
sometimes profitable to use diagrams of «,, and 1/@,, plotted against w; but they will 
not be required here. 


4. The Horizontal Line in Receptance Diagrams 


A node will be said to exist at a particular disc (the r'*, say) of a torsional system 
when that disc does not move during motion in a particular principal mode (the k™, 
say)*. In discussing this nodal condition, it is convenient to approach the state by 
a form of limiting process. 


Suppose that a torsional system is undergoing forced oscillation in all its 
principal modes according to equation (8), and the system is such that x, does vary 
during motion in the k“ mode. To discuss the condition when there is a node at x, 
(in this mode), a series of systems will be contemplated such that each differs slightly 
from its predecessor and such that the physical constants become so adjusted that 
the required nodal form is approached. This (ultimate) system is the one whose 
receptance functions will be examined. If exciting torques are applied at the 
co-ordinates x,, X,,...X, in turn and, in each case the response at x, is to become 
smaller in the k mode, it follows that the limiting process consists of letting .A,,, 
rAro,..-%Aen all tend to zero (see equation (8))f. 


Equation (6) may be written in the form 


cA re F 
_ (w?) + (7 — ©?) F (14) 
(w,? — g (w’) 
where F (w?) and g(w?) are functions of »?. On applying the limiting process to 
obtain a node at x,, so that ,A,, tends to zero, it is seen that 


(;? — ?) F (?) 


It follows, then, that all the functions f (w?) of equation (3) contain, as a factor, the 
k™ factor of A. These common factors do not affect the general shape of the 
receptance diagrams except where w=,. For this special value their ratio becomes 
indeterminate, as may be seen from equation (14), requiring a horizontal straight line 
to be drawn in all ,, and 1/2,, diagrams at the ordinate «,. 


(15) 


*This conforms with the usual usage of the word “ node.” Mathematically, this means that 
ox, /dp,=0, where p, is the k"® principal co-ordinate. The reader will perceive that this leads 
to an extension of the usual meaning of “ node” if x, is the r** generalised co-ordinate. 

+The shape of the k‘® principal mode of the limiting system is determined by the ratios between 
these constants that are approached as they all beeome small. 
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“rr 


(a) (b) 


Fig. 3. 


Yor, 


This result may be arrived at graphically. Taking the case of the direct recep- 
tance 2,,, the effect of the limiting process may be traced in Fig. 1. The k™ (dotted) 
resonance curve becomes very narrow at all points except those in the immediate 
neighbourhood of ,, since the numerator of the k™ partial fraction in the equation 
(6) for z,, becomes very small. In the limit, the curve is replaced by a straight 
horizontal line as shown in Fig. 3(a). The counterpart of this on the 1/¢,, diagram 
is shown in Fig. 3(5). 


The presence of the horizontal line in Fig. 3(5) means that the natural frequency 
of the k"* principal mode will not be altered if any sub-system is attached at the 
co-ordinate x,. 


Figure 4 shows a conveniently simple system which may be used to illustrate 
the result given in equation (15). If any of the receptances at x, are set up by the 
usual methods (e.g. see Ref. 3) it is found that the numerator has a factor in common 
with the denominator; this should be expected, since there is obviously a principal 
mode of oscillation in which the central mass remains stationary. It is found that 


—k(k+2A—2Mw?) (16) 
Mo? (k +224 —2 Mw?) (5k —2 Mw?) 

22°" Mo? (k+2A—2 Mw’) (5k —2 Mo?) 


a 


(17) 


The height, ,, of the straight line in Fig. 3(b) depends upon the system 
concerned. The line may occupy any position relative to the remaining branches 
and the rest of this paper deals with the two special cases that may exist. These 
are (I) where «, coincides with an anti-resonance frequency of the remaining 
branches and so lies on an asymptote and (II) where «, is a resonance frequency of 
the remaining branches so that the line intersects another branch at the w-axis. 
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5. Case I: w, an Anti-Resonance Frequency of the 
Remaining Branches 


If the point P in Fig. 3(a) lies on the vertical axis, no motion occurs at x, due 
to all the principal modes combined when the system is excited at x, with the 
frequency ,. This may be given an analytical interpretation as will now be shown. 


The series (6) for 2,, may now be written in the form 


= 


since it is required that the contributions of all the modes to the motion x, shall 
vanish as »—>o,. It follows that 


Opp = A (19) 


and, if the limiting process is now used to obtain a node at x, in the k™ principal 
mode (i.e. if ,A,,—> 0), then 


— 0) h(o?) 


—> 


That is, in the expression (3) for <,,, not only does f (w?) contain a factor that appears 
in A, but now that factor is repeated. 


2M 2M 
AVAVAVAVA 
k k 
Fig. 4. 


The system of Fig. 4 can be made to illustrate this result in an elementary way. 
By setting A=0, it is seen that (i) the cross-receptances of equation (16) retain the 
form of equation (15) whereas (ii) the direct receptance assumes the form of 
equation (20). 


The type of receptance function which is given by equation (20) can arise in a 
particular way. Suppose that a system may be divided into two sub-systems B and 
C such that x, is the only linking co-ordinate. It has been shown that the reciprocal 
receptance at x, for a system of this kind can be formed by adding those for each 
sub-system separately in accordance with equation (10). If this addition is performed 
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in a diagram like Fig. 2 when the two reciprocal receptances 1/8,, and I/7,, become 
infinite at nearly the same frequency, a branch is found which is nearly flat; as 
equality of these anti-resonance frequencies is approached, the horizontal line 
emerges at the ordinate w,. The composite system A will be capable of performing 
free oscillations with frequency ©,, the relative amplitudes of B and C determining 
the modal shape, since the forces at x, due to the two sub-systems must balance. 
It is also evident that if , is an anti-resonance frequency at x, for B and C, no 
motion will be caused at x, by a harmonic force or torque F,sinwt with this 
frequency. Since there is no motion in the k* mode due to the existence of the 
node, it may be asserted that the contribution of the remaining principal modes 
collectively is also zero. 


6. Case II: w, a Resonance Frequency of the Remaining 
Branches 


The case will now be examined where the intersection R of Fig. 3(b) lies on the 
vertical axis. This evidently implies that the natural frequency , (corresponding 
to the k™ principal mode in which a node exists at x,) coincides with an adjacent 
natural frequency ,. That is, A now has two equal roots. 


The direct receptance at x, may be written 


_ — 0?) — 0”) f, (©?) + (0? — ©”) + Apr — ©”)] f, (07) 
= 
If now o,, ‘ 
— [ f, ] 


The line appears in this expression together with the contents of the square brackets. 
The latter have the form of a direct receptance without a repeated factor and give 
rise to n—1 principal modes in the usual way. 


It is significant that, in the derivation of equation (21), use was not made of the 
limiting process ,A,,—>0 by means of which equation (15) was arrived at. The 
line which has been associated with a node has been arrived at apparently without 
the requirement that a node shall exist. Another new feature is that (21) is valid 
for all the co-ordinates; i.e. there is nothing in the derivation which refers to a 
particular co-ordinate x, The line appears in all the diagrams of reciprocal 
receptance and is always at the height of an intersection of the remaining branches 
with the w-axis. It is also seen that there is nothing which limits this argument to 
direct receptances; it is equally true for cross-receptances. 


These features can be illustrated by reference, once again, to the system of 
Fig. 4. If A=2k, the denominator of equations (16) and (17) is given the repeated 
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Fig. 5. 


root (Sk 2M’) and, when the receptance functions are found by the method of 
Ref. 3, the following results are reached : — 


whereas >= [ | | 


= ert —2/5 ] 
28° (6k—-2Mo*) L Mo? © 5k—2Mo? 


If attention is now confined to the partial fractions within the square brackets of 
equations (22) and (23) and principal modes are sought as in the derivation 
of equation (9), then those modes have the forms 


3:-2:-2 equation (22) (24) 
4:-—-1 by equation (23) 


Thus, the principal modes of the system of Fig. 4 (with A=2k) are no longer unique. 


In seeking an explanation of these new properties it is necessary to note the 
effect upon the principal modes of a repeated factor of A. If two principal modes 
have the same frequency, they lose their identities and one of them may be defined 
arbitrarily to be such as to provide a node at any selected co-ordinate. This follows 
because the system can then vibrate freely in any form which is a linear combination 
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of the two modes. Any such form can, alternatively, be taken as a principal mode 
and, in particular, one such form may be chosen which has a node at some specified 
co-ordinate. 


This may be illustrated simply by means of the system shown in Fig. 5, in 
which a particle performs small vibrations in the plane of the light springs. So long 
as the spring stiffnesses k and A differ, the system has natural frequencies given by 


(25) 


(for conformity with the notation, this implies that k< A). Each principal mode is 
obviously to be identified with motion parallel to the axis of one of the springs. 
However, if k is equal to A, the two natural frequencies are equal and the modes 
may be selected arbitrarily. Thus, they may be taken as distortions by which the 
particle moves along the lines at 45° to the axes of the springs; indeed, they can even 
be taken as motions by which the body describes circles in each direction about its 
mean position. 


FREQUENCY 
Fig. 6. 


This lack of uniqueness of the principal modes may be shown to exist by an 
analytical approach. It has been seen that counterparts of equation (21) may 
readily be found for all the ~,, and it was demonstrated in equation (24) that, after 
the extraction of the factors (w,?—*) from the top and bottom of this expression, 
n—1 principal modes may be found in the manner of equation (9). But it was 
found that the one having the frequency ©, is no longer unique since the ratios (9) 
now depend upon s. Let it be supposed that , and , are close, but distinct, roots 
of equation (4). If the excitation frequency is much closer to , than to any other 
natural frequency (apart from ,) as shown in Fig. 6, the receptances take the form 


rs iA rs 


“@,7-—0? 


(r=1,2,...m). (26) 


A limiting process is now contemplated in which the ratio a/b=C is preserved 
(where a=w,,— ©, b=,— ©, and C is some large fixed number) as a—> 0. Evidently 
the modal shape 


is obtained in this way. This has been found from a harmonic excitation applied at 
x, which must become vanishingly small as a—> 0 in order to retain a finite motion. 
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Had this excitation been applied at some other co-ordinate x, in this process a 
different shape would have been found because, in general, 


tAustiAus 


(28) 


Thus the modal shape depends upon the choice of the point of excitation x,. As has 
already been shown, this lack of uniqueness is to be expected on physical grounds. 


It has been shown that, if a node is connected ,to a fixed point through a 
spring, that is to say if stiffness is added at this co-ordinate, the frequency of the 
appropriate principal mode is left unaltered. Suppose that a system has a repeated 
factor (,’—w?) of A. The effect of adding stiffness at such a point is to render 
determinate the principal modes that were previously indeterminate; the one has a 
node at the point of attachment and has frequency ,, while the other has a higher 
frequency. As shown in Fig. 2, the effect of the spring is to shift the vertical axis 
of the (1/2,,, o) diagrams to the left, so that, here, the point R (Fig. 3(b)) is shifted 
thereby from the w-axis to a position such as is shown in Fig. 3(b). In this way the 
type of system is reduced to a less specialised form*. It is seen that, of the pair 
of frequencies which were originally equal, that one which is raised does not reach 
the adjacent anti-resonance frequency of the original system unless the anchoring 
spring is infinitely stiff. 


The converse of this process arises if one attempts to convert a non-specialised 
system to one having repeated factors of A. Suppose that one point only of the 
system is accessible and that this point is not a node for any principal mode. Then 
the diagram of the reciprocal of the direct receptance at this co-ordinate will not 
contain a horizontal line and the addition of mass or stiffness at this accessible point 
cannot create one (see Fig. 2). It may be deduced then, that the system cannot be 
modified in these ways so as to give its frequency equation repeated roots. On the 
other hand, if the accessible point is a node in one of the principal modes, it has 
been seen that the attachment of a spring at this point simply shifts the diagram 
of reciprocal receptance bodily to the right, so that an equal frequency system can 
be obtained in this way if the intersection R (Fig. 3) is to the left of the axis of w. 
If R is to the right of the axis, it may be brought into coincidence with it by the 
addition of mass at the accessible point; this second process distorts the receptance 
diagram as well as moving it. 


If two points are accessible, x, and x, say, then it may be possible to make 
x, a node by adding mass or stiffness at x,. After this, the system can be given 
equal frequencies by a modification at x, as before. Alternatively, a rigid, massless 
linkage may be made to join x, and x, such that some point within it is a node in 
one principal mode; now, by the previous argument, the equal-frequency system 
may be arrived at by a suitable modification at this point. 


*This may be illustrated by taking A=2k in the system of Fig. 4 (so that A has a repeated root) 
and then anchoring one of the end masses through a spring. 
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7. Conclusions 


The relation «,,=f (w*)/A gives rise to certain special cases. These obtain when 
the numerator and denominator have a common factor, and it has been shown that 
they imply the existence of a node at one of the co-ordinates. The node is associated 
with a horizontal line in receptance diagrams. The consequences of repetition of 
this common factor, either in f(*) or in A, have been examined. 


It has been assumed that, when a factor is repeated, it does not appear more 
than twice. This is not an essential requirement and its removal introduces no new 
concepts; it leads to the possibility of superimposing horizontal lines in receptance 
diagrams. Again, two different common factors of f(?) and A may be repeated 
and so on. 


The methods of this paper can be applied to generalised co-ordinates and, in 
particular, principal co-ordinates. In the latter event a node is said to exist at the 
r™ principal co-ordinate p, for all the remaining principal modes. This would imply 
that diagrams of the receptances at p,, p.,...p, Should each contain n—1 horizontal 
lines in addition to the single curve given by z= A/(w,?—’); the horizontal axis 
itself will be one of the straight lines if the system is not anchored. 
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On the Unsteady Motion of a Slender Body 
through a Compressible Fluid 


L. E. FRAENKEL 


(Department of Aeronautics, Imperial College)* 


SumMMaRY: Ward’s slender-body theory is extended to examine the unsteady 
motion of a slender body through a fluid at rest. The cross-sectional area 
and shape of the body, its forward velocity and its lateral motion are all 
arbitrary functions of time, but are subject to the restrictions of small 
disturbances. The length of the body is fixed. An approximate velocity 
potential is obtained by the joint use of two integral transforms, and the 
accuracy of this potential is discussed in some detail. General expressions 
for the aerodynamic forces acting on the body are derived in terms of a 
co-ordinate system which moves forward with the body in the mean direction 
of motion, but is fixed laterally. These expressions are then transformed, 
for the particular case of an oscillating rigid body moving forward with 
variable velocity, to expressions in terms of co-ordinates referred to the body 
axes. (The former co-ordinate system makes possible a fairly compact general 
treatment, whereas the latter is more convenient if the solution of a particular 
rigid-body problem is required.) Finally, arguments are advanced to justify 
the application of the slender-body theory to wings of slender plan form 
whose trailing edges are perpendicular to the direction of motion, and results 
are given for an oscillating wing moving forward with variable velocity. 


1. Introduction 


Since Ward’s slender-body theory"? has already been extended to problems of 
unsteady flow by Miles’ and Phythian™’, the appearance of a third paper on this 
subject requires some justification. Miles considered very general unsteady motions 
of a slender body, but confined his analysis to supersonic speeds. Phythian con- 
sidered both subsonic and supersonic speeds, but dealt only with particular types of 
motion. Miles worked principally with “wind axes,” whereas Phythian used “ bedy 
axes”; the transformation from one to the other is not always simple. 


In the present paper general motions at subsonic and supersonic speeds are 
considered and results referred to wind and body axes are correlated. The analysis 
is based on the joint use of two integral transforms and of the Stieltjes integral. This 
leads to what is believed to be a fresh discussion of the effect, on the accuracy of the 
slender-body solution, of discontinuities in the source and multi-source strengths 
and their derivatives. Correct general expressions for the drag are also derived for 
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*The work reported here was done while the author was in the Department of Aeronautics, 
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the first time: Miles’s result for drag is restricted to motions in which the body 
surface is fixed laterally*, while Phythian’s formulation of the momentum equation 
is incorrect for motions with rotary oscillations. 


It is extremely difficult to introduce this new matter without sketching in a 
suitable background, and accordingly the complete theory of the unsteady motion of 
slender bodies is briefly derived ab initio. 


NOTATIONT 
a,,%, Coefficients in expansion of slender-body potential 
b,,8, functions in slender-body potential 
b, b,,b, functions defining wing span (Section 8) 
B,B,,B, maximum values of b, b,, by 
c velocity of sound 
C,Cs body contours in plane s=constant, s=/ 
f= 
see equation (39) 
F _ see equation (53) 
g,y functions defining body surface 
h = 
I(t) distance travelled (see equation (9) ) 
J(s,t) see equation (34) 
k,«x doublet strength 
K, modified Bessel function of the second kind 
body length 
circumferential variable 
n normal outward from body in plane s=constant 
mn unit normal inward to volume of fluid 
p __ index of Laplace transform (Sections 3, 4, 8); static pressure (Section 5) 
q__ fluid velocity 


(r,6,s) cylindrical co-ordinates of frame moving longitudinally with body but 
fixed laterally 


*Professor Miles has explained in a private communication to the author that he calculated the 
drag under this restriction, but forgot to state it when summarising his results in Ref. 2. 


tCertain Greek symbols directly related to Latin ones, like z, (¢, 7)=a, (s, t), are defined in the 
first part of this list. 
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real part of 
R__ value of r on body surface 
S cross-sectional area of body 
S» body surface 
t time 
U(t) velocity of body in negative «-direction 
v see after equation (47) 
V_ volume of fluid 
Vn _—- see equation (59) 
w complex function such that Zw=9 
(x, y,s) Cartesian co-ordinates corresponding to (r, 9, s) 
z complex variable re 
a,B8 body incidence in xs, ys-plane 
€ maximum value of R 
p density 
o longitudinal co-ordinate from fixed origin, 
r time 
¢,¢  Sslender-body velocity potential 
errors of see Appendix 
® exact velocity potential 
w index of Fourier transform 
angular velocity of rotating frame 
( ), conditions at infinity 
( ), quantities referred to body axes 


2. Formulation of the Problem 


Let (x, y, «) be a right-handed system of Cartesian co-ordinates fixed in a fluid 
at rest at infinity, let (r, 0,0) be the corresponding system of cylindrical polar co- 
ordinates, and let z denote the time. A slender body, whose nose and base rest at 
o=0 and o=l, respectively, until time s=0, begins to move at this instant in the 
direction of the negative o-axis with velocity U (7), small lateral departures from 
the o-axis being permitted. 


Let « denote the upper limit of r for points on the body; we specify that the 
length of the body be O(1) and that the Mach number of the forward motion and all 
time rates of change be bounded (the forward and lateral velocities of the body may, 
however, be discontinuous): the case «—>0 is then examined. 
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The words “cross section” and “body contour” are used to describe sections 
and the body boundary in planes normal to the longitudinal body axis, and “ axial 
section ” and “ profile ” are used to describe sections and the body boundary in planes 
through this axis. The body is subject to the restrictions specified by Ward, except 
that its cross-sectional area and shape may vary with time and that, initially, bodies 
are considered which fail to satisfy, in regions whose total length is at most O(c), 
the requirement of small profile slope and curvature: that is, the effects of rounded 
and open noses and tails and of discontinuities in profile slope are considered 
initially. Rounded noses and tails are assumed to be such that the cross-sectional 
area and its derivatives are everywhere O (ce), (this requirement is satisfied, for 
example, by a slender ellipsoid). 


Consider the irrotational flow of a perfect, inviscid fluid. Let ®(x, y,o,7) be 
the velocity potential, c the speed of sound at any point and time, c, the speed of 
sound in the undisturbed fluid, and let T° be the ratio of the specific heats. The 
governing differential equation is 


a ] 


where e=c,'— +vay 


(1) 


If © and its derivatives are small, equation (1) may be linearised to give the wave 
equation 


1 


3. The Slender-Body Solution 


We introduce Fourier and Laplace transforms: — 
F(x, y.0, t)= [em F (x,y, a,t)do, (3) 


which we also write F(x, y,w,7) =F (x. y,¢,7), 
and F(x, y,, p)= | per (x.y. 0,7) dr, 


or F(x, y,0,p)# F (x, 7). 


To obtain the subsidiary equation of (2), we multiply it by the kernels of the 
integrals and integrate. If there are points (x;,y.,0.,7,) Where the solution will 
prove singular and (2) itself is not satisfied, the subsidiary equation cannot be 
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justified on lines x=x,, y=y,; but if elsewhere (i) ® and its first partial derivatives 
are continuous with respect to space and time, and (ii) ® and its first and second 
partial derivatives —> 0 as «—> oO in such a manner that integrals of the type (3) 
converge, it is found, upon integrating by parts, that the subsidiary equation is 


eb ah 
—h? (w, p)®, 


where h? (w, p) =? + p?/c,?. 


If the first partial derivatives of ® are discontinuous across spherical surfaces of 
the form 

Co” (7-7, P=(X- + Yo)” +(¢-<¢,)’, 
or across surfaces which consist of an envelope of such spheres, (® itself being 
continuous), & still satisfies (5); for it may be shown by the method of Ref. 4 that 


the sum of the additional terms which each discontinuity contributes to the 
subsidiary equation is zero. 


A general solution of (5), representing disturbances on the c-axis, is 
(r,0, 0, p)= A, p) Ky (hr) + K, A, (o, p)cosn6+ B,(, p)sinnd], (©) 


where the K,, are modified Bessel functions of the second kind’, and the A,, B, are 
arbitrary functions. It is assumed for the moment that » and p may be treated as 
O(1); then since r is O(¢) on and near the body only the dominant term in the 
series expansion of each Bessel function is retained, the apparent error being a 
factor 1+O/(«?). There results the approximate potential 


where z=re®, &=-A,, B,=&(E+logh/2), 
(8) 
&,=2"-!(n—1)! h-* (A, +iB,), 


and E is Euler’s constant. Like Ward’s solution for the steady supersonic case, 
this potential is harmonic in (r, 4); the coefficients <, are determined by finding the 
harmonic function which satisfies the boundary condition and expanding it in series. 
The series converges outside a circular cylinder enclosing the body. 


We now introduce variables s, t such that s is measured along the c-axis from 
the plane of the nose of the body. We write 


U (7) d7’ =I (r)=1 (9) 


s=c+I(r), ft=r, . (10) 
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and let 0/0t imply r,9, s held constant, whereas 0/07 implies r,6,o0 held constant. 
Equation (7) may now be written 


¢(r, 9,5, F {a, (s, t) log z+ b, (s, 2 


where a, (s, t)=2,(o,7), ete. 


4. The Evaluation of a. and b, 
If the surface of the body is defined by 


7)=g(x, 5, t)=0, 


the boundary condition is Vy. Ve= on the body, 


This boundary condition, together with our restrictions on the shape of the 
body and a knowledge of two-dimensional incompressible flows, ensures that 09/0x 
and 0¢/0y are O(«) everywhere except in the immediate neighbourhood of a 
rounded nose or tail, (where 0g/0s is larger than O(¢)). Hence the a, are O(e"*’), 
except in these regions; ¢, 09/0s, 0¢/0t, etc., are O (e* log ¢); and 079/0x*, 0°o/dy’, 
etc., are O(1). Accordingly, equation (12) becomes 


ax Ox * ay ay LUO ag 
Op on On 


where 7 is the normal outward from the body contour in a plane s=constant, and, 
by a special notation, 


on dg/ on dg/ Ot 
ato {(0g/ Ox)? + (0g/0y)*} 


(on/ds and On/Ot are defined as positive when the boundary moves outwards with 
increasing s or 


(14) 


By integrating 0¢/0n around the body contour in a plane s=constant, using 
first (11) and then (13), we find that 


a, (s, 


sin, 


where S(s, t) is the cross-sectional area of the body. In the regions s<.0, s>/, or 
t <0, a, and all the a, vanish. 
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Bo (, p) is now inverted by recalling that 
B)+4,logr=4, (E+1og =) =lim,_,), — (Ar). (16) 


K, (hr) may readily be inverted by using the identity’? 


and inverting with respect to w first. It is found that 


] 24 


K, (hr) = 


where H is the Heaviside unit function. This is the we.l-known acoustic source 
potential. The double Faltung integral appropriate to our integral transforms is 


Foo, p)G(o, do’ | 18) 


r= 


where the integral with respect to 7’ is a Stieltjes integral; so that in the present case 


(hr) 4, = ax’ | da, :’) 
=0- 


By means of the rule for integration by parts based on the Stieltjes integral, this 
becomes 


~K,(hr)é, = | sinh- 2 


o=-c 


where we now have a Stieltjes integral with respect to o’. The relation 
a, (A, B)=a,(A+I1(B), B) is used, and the variable of integration is changed from 
a’ to 7’, where 
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+ | sinh Vice da, ( + 


c 


r 


=— 5 | cosh [ da (s— Vie 


Hence B, +4, log r=lim —K, (hr) 4, 


t 


1 2c, 4 é J 


+ da, (s+ +10, (23) 


t 


and b, (s, t)= — a, (s, t) log 2c, — | log (t—1’) [ da, (s—c, (0), r) 


+ da, (s+e, r)] » 


Some physical significance may be attached to the foregoing expressions by con- 
sidering an (r, «, 7)-space (Fig. 1). A disturbance at a point (0, 0’, 7’), which lies between 
the trajectories of the nose and tail of the body, affects only points (r, «, 7) on its “after- 
cone,” c,(r-7’)= +r°*}. Conversely, in an (r’,o’,7’)-space the flow at a 
point (r,o,r) is affected only by disturbances (0,0’,7’) lying on the hyperbola 
Cy (7’ —7)= — which is formed by the intersection of the “ forecone ” 
of (r,s, t) with the plane ”°=0. If all the aftercones of disturbances were drawn in 
Fig. 1, a section by a plane s=constant would indicate the disturbed field at that 
instant in a physical, ro-plane. The slender-body approximation assumes r— 0+. 
and replaces the cones by planes tangent to them and parallel to the r-axis. 


If U(t)=c,, the function 5, (s,t) tends logarithmically to infinity as t—> 0. 
that is, as the steady state is approached. The other terms in the slender-bod) 
solution remain bounded and O (<? log ¢). 


66 The Aeronautical Quarter! 


Hence 
r’ = 
— 0 
5. 
bx 
fu 
r 0 )+ D 
st 
af 
in 
O 
t'=0- a 
ve 
: 
de 
ve 
to 
di 
Fe 


22) 


23) 


ter! 


UNSTEADY MOTION OF A SLENDER BODY 


NOSE TRAJECTORY 


| 


TAIL TRAJECTORY 


Fig. 1. 


5. The Drag Force 


The accuracy of the slender-body solution is discussed in the Appendix. As a 
result of the conclusions drawn there, the following restriction is now applied. The 
body must be such that the source and multi-source strengths a, (s,t) are continuous 
functions of s everywhere, except at the base of the body, s=/, in supersonic flow. 
Discontinuities with respect to time are permitted in the source and multi-source 
strengths, but it must be understood that the results do not provide an acceptable 
approximation at instants immediately following such discontinuities. Discon- 
tinuities are also permitted in the 0a,,/0s; it is shown-in the Appendix that these 
increase the errors of the integrated forces given in what follows by a factor of 
O (log £) at most. 


To calculate the aerodynamic forces, consider the flow of momentum through 
a surface consisting of the body surface §, and a circular cylinder moving with 
velocity U (t) along the -"-axis in the negative direction (Fig. 2). The cylinder has a 
small radius; its curved surface is denoted by S,, and its plane ends at s=0 and 
s=l are denoted by S, and S, respectively. The contour of the body at the base is 
denoted by Cy. For such a surface S, containing a volume V, (the boundary varying 
with time), the momentum equation may be written 

d 


{ (p-p,)nds= J ota -w).njqdS,. 
S(t) V(t) S(t) 

where p is the local pressure, p the local density, q the local fluid velocity, w the 

velocity of any point on the control surface, and n is the unit normal drawn inwards 

to the volume of fluid. Hence, taking the s-component of (25), and ignoring the base 

drag, we obtain 


5,-8, P(t) 5, 
+O(e*log?<). . (26) 
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( 
u(t) Sb 


Fig. 2. 


From the energy equation we have 


(27) 


1+0 (+ log). 


0 


Hence, noting that 09/dr=09/rd0=0 in S,, we obtain 


53 


- 2 4154. O(c log?®). . (28) 


or Os 
Sy 
4 (2% gy, 
(t t 


and by Gauss’s theorem, writing k for the unit vector in the s-direction, we have 


ot 
V(t) S(t) S,-S, 3 Sy 
Hence di dvV= dS + th ht ds. (9) 
(t) S,-S 
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The following result is taken from Ward’s analysis (Ref. 1, equations (34) to (36) ) 


1 1 2% _ ab, 5 2 

24 + r 06 "Jas or 95 S=* ds— 5 
0 


(30) 


2 


where m denotes the circumferential variable. By (29) and (30), (28) now becomes 


0a, 0b, On On 
0 


CR Sp 


(31) 


This result is valid for both subsonic and supersonic speeds, but the body is permitted 
to have a base only for supersonic speeds. 


Finally an identity is derived which will prove useful later. By Gauss’s theorem 
~ dt = as | 


V(t) S,-S, Sp 


Comparing (29) and (32) gives 


an _ 29 _ an 


5, Sy 


6. The Lateral Forces 


The lateral forces may be calculated by following Ward’s analysis and by dealing 
with the unsteady flow terms in a manner analogous to that of the previous section. 


The function J(s,j=-i pdz. ‘ . (34) 


is introduced, where C denotes the body contour in a plane s=constant and z, is the 
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complex co-ordinate of the centroid of area of the body cross section. The complex 
lateral force is then given by 


l 


0 


0 


0 0 


The complex moment is (M, and M, denote the moments about the x- and y-axes 
in the positive sense for a right-handed system) 


0 


7. The Oscillating Rigid Body 


Let (x,.y,,5,,t,) be variables referred to the axes of a rigid body which is 
moving with velocity U (t,) in the direction of the negative c-axis and is oscillating 
with angular velocity —.’(t,) about an axis x,—0, s,;=s,.(¢,), and with angular 
velocity about an axis y,=0, s,=s,»(t,). The angles and are taken as 
positive when the body nose points in the positive x and y directions, respectively. 
Since it is easier to solve the two-dimensional harmonic problem in terms of these 
new co-ordinates, we wish to transform our expressions for drag and lift into 
expressions referred to the new system. For points near the body the new frame is 
related to the previous one by the transformation 


x, =x+ +O +0 
(38) 
s,=s+0O(e?), 

0/0t, implies x,, y,,s, held constant. The potential may be written 


where, by expanding log z, z-', etc., in terms of z,, and writing f,+if,=f, we find 
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Qi (s,, t,)=a, (s, t), b,, (s, 
(41) 


4,, (S,,t,)= —a, (s, t) f (s, (s, 0), etc. 


Let the surface of the body be defined by g,(x,,y,.5,)=0. Then by (14) 
and (38) 


os as, Ox, Os oy, + 
or Ox, Ot dy, Ot dy, 


In the expression for drag we have, using (33), 


(42) 


b 


where (cf. equation (14)) —0n,/0s, is the approximate s,-component of the unit 
normal to the body surface, and is independent of the time. Recalling the definitions 
of the functions J, f, and f, (equations (34) and (39)), and writing «+i8=n, we 
obtain for the expression (43) 


Ss 0 


(t,) )- iB’ (t,) ) (44) 


and 


Po 


0 


= 1) 500) +18 +O log? <). . (45) 
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The expressions for lift and moment are unchanged; (s,t) may be replaced by 
(s,,t,). It remains to express J in terms of a,, and z,,. By equation (35), we have 
for the present problem, using equations (41), (15) and (38), 

J (s,,t,)=2 7a, (s, t)+ \[sozo] 


U (t,) S’(s,) 
] 


+ + | [ | 


=2" [ a, (s,,.t,) +f 


1 


The problem treated in this section is that considered by Phythian®. Equation 
(45) corrects Phythian’s expression for drag (i.e. his equation (72)), but does not 
reveal the cause of his error. Phythian used a control surface of constant shape, fixed 
to the body and rotating with it, and variables corresponding to our (x,, y,,5;,1f,). 
He took components of momentum in the instantaneous directions of Ox,, Oy,,Os,. 
In such an approach the volume integral in the momentum equation (25) becomes 


where @ is the angular velocity of the rotating frame. Phythian omitted the vector 
product term; if this is added to his analysis the result given here is obtained. In 
principle the extra term also contributes to the lift, but this contribution proves to be 
of negligible order, so that Phythian’s expression for lift is correct. 


8. Laminar Wings of Slender Plan Form 


We approach the problem of laminar wings of slender plan form by consider- 
ing first the ordinary linearised solution for a wing having a span of O(1) and 
moving forward approximately in the yo-plane, its displacements and velocities due 
to pitching, plunging and rolling being O(¢). The wing and its wake are represented 
by a distribution of doublets; from the viewpoint of the ordinary theory it is 
immaterial whether these are placed in the yo-plane or in the instantaneous plane of 
the wing but, since slender-body theory requires the boundary condition to be 
satisfied on the body surface, the latter course is chosen. We denote variables 
referred to the wing axes by (x,, y,,5,,t,) and define o,=s,—J(t,), 7,=t,. (0/07, 
implies x,,y,,07, held constant.) The starboard and port edges of the wing are 
x,=0, y,=—b,(s,) and y,=b, (s,), respectively, the maximum values of b, and b, 
being denoted by B, and B,. To keep the slender-wing doublet strength k and the 
coefficient a,, continuous, we specify that the wing be pointed, b(0)=0. (See 
equation (54) ). 


72 The Aeronautical Quarterly 


| 
i 
J 
O 
i 
F 


UNSTEADY MOTION OF A SLENDER BODY 


0 0 


where v is the velocity of the origin of the (x,, y,,7,) frame, && is the angular velocity 
of this frame, and r, is a position vector with components x,,y,,7,. v and @& are 
O(«), and the linearised equation may therefore be written 


= 1400), (48) 
where the neglected terms are of the same type as those neglected in equation (2). 
Hence, if », and p, are Fourier and Laplace operators corresponding to «, and 7,, 
and h,?=w,?+p,?/c,?, our solution is 


where «(y,,0,,7,)=k(y,,5,,¢,) is the doublet strength. Near wings of slender 
plan form, and in the absence of discontinuities in «, the argument of the Bessel 
function may be treated as O(«); accordingly only the first term in the series 
expansion of K, is retained. The validity of this step is discussed in what follows. 
There results the approximate potential 


b 


x 
dy 
Ok (y,’,5,.t,) , 
= ( ay,’ tan dy,’. ($2) 


This potential is again harmonic in (x,,y,): the solution of the integral equation 
which determines k for a given normal velocity on the wing is well known. However, 
if there are lines s,=constant, x,=0, which are occupied partly by the wing and 
partly by the wake, substantial difficulties are introduced. These are avoided by 
considering only wings whose trailing edges are straight and normal to the direction 
of motion. 


It may be noted in passing that rolling velocities and displacements of O (1) 
can now be admitted for a wing of slender plan form, since these only cause velocities 
of O(<) normal to the wing surface, and since the solution on the wing is now 
independent of conditions in the wake. 
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To illustrate the remarks that follow, consider a wing of slender symmetrical 
plan form (b,=5,=5), oscillating in pitch like the body of the previous section, with 
B(t,)=0. The boundary condition is 


(t,) a(t,)— 2’ (t,)[s, — Si (t,)) = —F(s,,t,), say. (53) 


Hence [ (s,)} -2) |, 


(54) 


FF (5.11) 


It is apparent that 0k/0s, is not only discontinuous but also infinite at the lead- 
ing edge (Ward’s restriction on the curvature of the body contour has been violated); 
hence the approximation to (49) is likely to be less accurate than for bodies. The 
singularity is, however, of the same order as that found in more exact solutions, 
and the integrated errors are generally assumed to be acceptable. Ward") has 
discussed this point fully. 


At subsonic speeds the Kutta-Joukowski condition should be satisfied at the 
trailing edge, and it is evident that the present solution fails to meet this require- 
ment. This may be attributed to a finite discontinuity in 0k/0s, and da,,/0s, across 
the trailing edge (k and a,, remaining continuous). But it has been seen that this 
type of discontinuity is admissible on a slender body in spite of the large localised 
errors it introduces, and there is no reason why the integrated error should be 
greater in the present case. 


The expressions for the aerodynamic forces on the oscillating wing are readily 
obtained. We have 


a,,=6,,=0, &,). 


F? (1, t,) B?. (55) 


B 
tm, = 2 | VB dy,= 


Hence from equations (45), (36) and (37), 


x B?U [ | +¥, | +V, 2", a) 


0 


| -V, [ +V, [ | + 
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where | s,"b*(s,)ds,, (m=0,1,2).. . . . (59) 
0 


Because in the present problem we have a,=0, a,=a,,, and an anti-symmetric 
potential, these results would also have been obtained.by satisfying the boundary 
condition in the fixed yo-plane. Such a procedure would yield correct pressure 
differences across the wing, but incorrect pressures (see equations (27) and (47)), 
and it would be inconsistent with the procedure which must be used for bodies 


with thickness. 


The results of this section have been largely given by Miles, who used a 
somewhat different approach and who also questioned the validity of the method 
for subsonic speeds. 
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Appendix 
THE ACCURACY OF THE SLENDER-BoDy SOLUTION 


(i) The Accuracy in the Absence of Discontinuities 

If the a, (s, t) and their first partial derivatives are continuous functions of s and t 
everywhere, the error of the slender- body potential is uniform, for values of r on and 
near the body, and can be found by using Ward’s technique"). In terms of z=x+iy, 
Z=x-—iy, o and 7, we have 


and the exact equation (1) may be written 


To obtain a second approximation, put 
P=o+¢9 +[w (z)+ w w(z J+ $1 


where w is the complex function of which the slender-body potential is the real part, and 
¢y Tepresents the higher order terms to be added. It is assumed that the contribution 
of y, to the right hand side of (61) is negligible compared with the contribution of ¢. 
The right hand side of equation (61) is then known, and an appropriate form of 
solution is ¢;=9,+¢,, Where ¢, is a particular integral of the equation 


and ¢, is a solution of the linearised equation, having the appropriate properties at 
infinity and satisfying the boundary condition 


Ward’s method may be used to show that ¢,, 09,/0c, 09,/0r, etc., are O (s* log’ ¢), 
while 0y,/0z, etc., are log It follows that these are also the orders 
of ¢, and its derivatives. It is worth noting that if terms not appearing in the linearised 
equation had been omitted from the foregoing analysis, g, and gy, would have been 
found to be O («* log <), which is very nearly correct. 


(ii) The Effect of Discontinuities in the a, (s, t) 


To examine the localised errors introduced by discontinuities in the a, (s,f) and 
their derivatives, the accuracy of the slender-body potential is considered only as a 
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solution of the linearised equation. In what follows the slender-body solution is 
denoted by ¢ (x, y, o, 7) or @ (x, y, s, 1); the exact solution of the linearised equation with 
the source and multi-source strengths of the slender-body approximation, (equation (6) 
with the A, and B,, defined by equation (8) ), is denoted by ¢+4,; and the exact solution 
of the linearised equation which satisfies the boundary condition (13) on the body surface 
is denoted by ¢+4,+ 9p. 


If a, (s, t) is discontinuous at s=A for all t, (e.g. if the body has a rounded or open 
nose or tail, or a discontinuity in profile slope), then in 9+ 9, there are two terms, (by 
equation (22)), 


where the ¢,’ are defined by 


SEV 7. 65) 


0 


In the slender-body approximation ¢ the corresponding terms are, (by equation (23) ), 
2c, (t-t 
(A+, - a, ts,’ | log fi) 


Now it can be shown that in the region |s- A| > O(¢), c, | ts’—ts,’| is O{e?/(s—A)}, 
provided that U (t’) is not near sonic speed in the interval between ¢,’ and ¢;,’. (If 
U (¢’)=constant this restriction is not necessary). Hence in the region |s- A| > O(c), 
is 


e4 


¢» is closely related to 9,. By definition ¢, is an exact solution of the linearised 
equation such that on the body 0,/0n= — 09,/0n and such that ¢, and its derivatives 
vanish as r,s—> 00. In view of these conditions and of knowledge of the order of 9,, 
it seems safe to assume that in the range |s- A|>O(e), ¢, is at most 


€ 


This assumption is confirmed by such more exact solutions as are available; these being, 
for the case of a rounded nose at subsonic speed, the classical solution for an ellipsoid 
of revolution, and for the case of a discontinuity in profile slope at supersonic speed, 
Ward’s quasi-cylinder solution. 


If a, (s, t) is discontinuous at t=B for all s on the body, (e.g. if the body changes its 
forward velocity discontinuously or expands discontinuously), the effect is similar. For 
c, (t - B)>O (¢), and for all trajectories, ¢, and ¢, are of the order given by equations 
(68) and (69) with |s—A| replaced by c, (t— B). 
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These observations must be modified when we consider the potential at points (r,r,7) 
or (r,s, t) lying between the after-cone of the point of intersection of two discontinuity 
lines s’= A, t/=B, and the vertical tangent plane of this cone (Fig. 3). The contribution 


Fig. 3. 


of the s’ =A discontinuity to ¢+¢,, which springs from the point S, then differs from the 
contribution to ¢, which springs from S,, by an amount greater than that given, for it is 
no longer true that 


2 
a, (At, ts,)=a, (A+,ts)4+0 
s—A 


The same applies to the contributions of the “=B discontinuity. However, there exist 


instead relations like 
2 


so that the combined contribution to ¢+4, of the points § and T is still equal to the 
contribution to » of the points S, and 7, within the errors given. A similar argument 
applies when the body changes its velocity discontinuously at  =B from a substantially 
subsonic to a substantially supersonic value. 


The effect of discontinuities in the multi-source strengths a, (n > 1) can be studied 
similarly. It is found that the errors are of the same order as those already quoted, 
except that no log terms appear. 


The expressions (66) may be written 


Utc, Utc, (Utc,)r 


—T(ts,’) 


t—ts,’ 


(70) 


where U= 


It is evident that such terms in 9, if used in the region s- A < O(e) where they are 
invalid, lead to a non-integrable singularity in pressure at s=A, (of O {(s—A)~-'}), and 
therefore to infinite drag. Furthermore, it is known that in a more exact solution the 
region of a body close to a rounded nose or to a discontinuity in profile slope generally 
makes contributions to the aerodynamic forces which are of the order of these forces 
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themselves. For supersonic speeds it is possible partially to overcome this difficulty’ *, 
but in the present paper the forces are calculated only for bodies satisfying the 
restriction mentioned in Section 5. 


(iii) The Effect of Discontinuities in the 0a,/0s and 0a,/ Ot 


That part of the source potential -K,(hr)&@, which is not the result of 
discontinuities in 2, may be written , 


oo 
1 


«x 
1 o’—o 1 )] 
OF (o, oF 
where the notation OF (¢,7) 
Oo 


has been introduced. 


The expression (71) represents a sum of Riemann integrals taken between the points 
of discontinuity of the integrand. We can again apply Stieltjes integration by parts to 
(71), and transform to the variables (s, t), in order to examine the contributions of dis- 
continuities in da,/ds and da,/dt. A similar procedure can be applied to the multi- 
source potentials. 


It is found that if da,,/ds is discontinuous at s=A for all ¢, (e.g. if the body has a 
discontinuity in profile curvature at s=A), ¢, is 


of 7 for |s—A|>O(e) (72) 


s— 


and o| for |s-A|<O(e), . (73) 


log 


where the log term appears only for n=0. It would be sufficient to write (73) as O (¢*); 
the extended form has been written to indicate the behaviour near s= A of ¢+4,, which 
contributes the first term, and of ¢, which contributes the second. It is evident that 9 
and ¢, are continuous at s= A, but that errors in pressure can be O (<* log|s—A|) and 
logarithmically infinite. For the reasons already stated, it is assumed that (72) and (73) 
also give the orders of 9. 


If da,,/0t is discontinuous at t=B for all s on the body, the orders of the error are 
given by (72) and (73) with |s—A| replaced by c, (t—B). 


On the basis of these results it is shown in what follows that, in spite of the 
substantial localised errors in pressure, discontinuities in the da,/0s only increase the 
error of the integrated expressions for the aerodynamic forces by a factor of O (log «), 
provided that the discontinuities are separated by distances of O (1). 


The foregoing assumption about the order of , is confirmed by the fact that the 
slender-body approximation for the steady subsonic flow past a body pointed at both 
ends (the body must be pointed if the a, are to be continuous functions of s) yields the 
correct values of zero lift and drag, even though the 0a,/ds be discontinuous. For super- 
sonic speeds the assumption is again confirmed by the quasi-cylinder solution. 
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(iv) The Effect on the Aerodynamic Forces of a Discontinuity in the 0a,/0s 


Let the 0a,/0s be discontinuous at s=0 and let the body extend from s=-—!/, 
to l,. Let ¢; denote the error of the slender-body solution. Then it follows from the 
combined results given earlier that there exist positive finite numbers C,,C,,C, such that 


for |s| «© « (4) 
ef 
and for |s| >Cye, $1<.C, Ts] log? Ts] (75) 
The effect of this error on the drag force is given by 


3(5%) 
Os” “Ot or or or 


R? 06 00 + 


where r=R (6,5, t) denotes the body surface. Now 


and ¢, is a continuous function of s. Hence, we have from the 0¢,/ds term in 
equation (76), 


Cre 
Ss 
=O {« log? «+ | e'ds+ | log’ <as } 
0 
=O(e® loge «). . ‘ . (78) 
The remainder of (76) leads to terms of this order or smaller. The effect upon 


the lateral forces is similar: for the force X, we merely replace R (0R/ds) in (76) by 
(09R/06) sin 6+ R cos @. 
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